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SUMMARY 


The  impingem  nt  of  lie  I  ^yt^rgolic  propellants  following  mixing  and 
separation  mechanisms  is  characterized  by  chemical  reaction  time  and  residence 
time  criteria.  Experimental  results,  obtained  by  hydrazine  impinging  with 
nitrogen  tetroxlde,  agree  very  wall  with  the  limit  defined  by  assuming  that 
separation  occurs  when  residence  time  is  greater  than  reaction  time.  This 
criteria  shows  a  strong  influence  of  propellant  temperature  because  reaction 
time  is  exponentially  dependent  upon  temperature  in  this  ignition  region  (the 
existence  of  an  ambient  tsmperature  ignition  region  is  the  essence  of  the 
meaning  of  hypergollcity) .  Further  experimental  results  indicate  that  chamber 
pressure  had  only  a  second  order  effect  and  that  monomechyl  hydrazine  and 
TTDM*v/hydrar*n<*  (CO/SC  blend)  impinging  v.Mth  nitrogen  tetroxlde  separate  at 
higher  propellant  temperatures  than  those  obtained  for  the  hydra zine~nltrogen 
tetroxlde  combination.  Stream  separation  always  occurred  when  impinging 
hydrazine  with  compound  A  (C*F»)  in  the  experimental  fuel  and  oxidizer 
temperature  range  (hydrazine  -40  to  60°F  with  compound  A  -  0  to  30°F). 

A  second  work  area  of  the  project  concerned  the  extension  of  the 
Steady-State  Combustion  Computer  Model  undert  ;en  as  a  joint  technical 
effort  between  the  Air  Force  Rocket  Propulsion  Laboratory  and  Dynamic 
Science.  The  computer  program  was  generalised  to  allow  increased 
versatility  in  running  various  propellant  syetoms.  The  numerical  procedures 
end  organization  of  the  program  were  substantially  Improved.  Propellant 
property  values  were  gathered  for  several  sj  items  of  current  interest.  A 
computer  program  listing  is  presented. 
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PART  I.  REACTIVE  STREAM  IMPINGEMENT 


1.  INTRODUCTION 

The  work  reported  here  Involves  an  experimental  determination  of  the 
physical  and  design  parameters  which  control  the  mechanisms  of  hypergollc 
propellant  stream  Impingement.  The  need  and  usefulness  for  this  work  Is 
discussed  and  a  physical  picture  of  hypergollc  Impingement  was  developed. 

This  picture  led  to  a  model  of  hypergollc  stream  mixing  reaction  and  blow- 
apart  criteria  and  the  criteria  were  used  to  define  an  experimental  technique 
Involving  photographic  observation  of  the  Impingement  plane.  Mix/separate 
criteria  were  measured  for  nitrogen  tetroxide  impingement  with  hydrazine,  MMH, 
and  50-50  (hydrazlne-UDMH) .  Compound  A  (Cl  Tb)  was  also  Investigated  and 
photographically  observed  The  results  of  :he  experimental  criteria  were  applied 
to  demonstrate  their  effect  on  rocket  engii  a  design  and  were  compared  to  large 
scale  engine  performance  and  injection  mixing  explosion  measurements. 

a  *  Propellant  Injector  Design 

Injection  spray  patterns  control  O/F  mixing  and  mass  distribution 
within  a  rocket  engine  combustion  chamber  and,  consequently,  these  spray 
patterns  exert  a  major  design  influence  over  chamber  performance,  thrurt  chamber 
liner  compatibility  and  combustion  instability.  Injector  design  was  largely  a 
matter  of  intuition  and  trial  and  error  testing  until  the  early  work  of  Rupe,(Ref.l) 
demonstrated  the  possibility  of  defining  the  relationship  of  design  parameters  to 
the  uniformity  of  stream  mixing.  The  mixing  uniformity  criteria  were  developed 
by  Ruptr  for  nonreactlve  impingement  and  involved  two  phases  of  injector  design: 
first,  the  spray  pattern  emanating  from  the  individual  set  of  impinging  elements 
must  be  defined  in  terms  of  the  design  parameters  which  affect  impingement 
characteristics  and,  secondly,  all  the  spray  patterns  together  must  be  layed-out 
to  give  the  mos*  uniform  mass  flux  within  the  constraint  of  the  design  manifolding. 
Design  of  unlike  hypergollc  Injectors  has  remained  a  matter  of  trial  and  error, 
however,  because  the  '  st  phase  determination,  that  of  the  spray  pattern 
emanating  fiom  the  Individual  set  of  unlike  hypergollc  impinging  elements,  has 
not  been  defined  in  terms  of  the  design  operating  parameters  which  affect 
impingement  characteristics . 


b.  Hypergollc  Impingement 


The  question  which  must  be  answered  before  trial  and  error  testing 
of  hypergollc  injectors  will  lead  to  any  fruitful  Lntuitlon  Is:  Do  hypergollc 
streams  mix  nonreactlvel/  or  do  they  always  separate  (t.e.  ,  Plow  apart),  or 
perhaps  for  different  critical  stream  parameters  ir  there  mix/separate  criteria  ? 
Previous  work  has  Indicated  that  hypergollc  streams  sometimes  mix  and  sometimes 
separate.  It  was  the  purpose  of  this  work,  therefore,  to  experimentally  define 
which  parameters  are  most  influential  in  defining  the  mix/separate  phenomena 
of  several  hypergollc  combinations . 

In  any  impingement  and  atomization  scheme  the  mixing  uniformity  would 
appear  to  be  important  and  it  can  be  related  to  a  ratio  such  as  the  streams' 


momentum  or  stagnation  pressures  which 
orifice  diameter,  D.  Now,  for  extension 
of  hypergc! icity  must  be  considered  as  it 
influenced  cy  temperature,?)  and  re s1- den 


involve  the  density,  o  ;  velocity,  v;  and 
to  hypergollc  streams,  the  basic  nature 
involves  chemical  kinetics (exponential! 
re  time  within  a  heat  balance  control 


voiume.  .Another  system  parameter  which  effects  chemical  kinetics  is  chamber 
pies  Sure ,  while  orifice  hydraulics  (L/D)  always  affect  turbulence  l^vcl  and 
mixing.  Having  listed  these  parameters,  it  is  easily  seen  why  trial  and 
error  engine  testing  of  hypergolic  impingement  efficiency  and  popping  has 
nor.  been  fruitful:  the  chemical  kinetics  of  l.ypergolicity  are  exponentially 
temperature  dependent  while  in  full  scale  engine  test  programs  this  parameter 
is  seldom  considered;  it  is  often  not  controlled  from  run-to-run;  and  is  also 
often  .not  even  measured  or  available  in  data  reduction. 

c.  Previous  Work 

Bari/  investigations  in  this  field  have  tended  to  be  phenomenological! 
oriented,  thus,  the  work  reported  by  Elverem  and  Staudnammer  (Ref.  2),  and 
Johnson  (Ref.  3)  noted  that  streams  do  separate  and  they  continued  to  show  the 
magnitude  of  the  effect.  The  later  work  of  Burrows  (Ref.  4)  gave  indication 
that  hypergolic  streams  do  not  always  separate  and  measured  specific  conditions 
necessary  for  separation  to  occur.  More  extensive  work  has  been  carried  out 
oy  lawver  and  Breen  (Ret.  5)  to  define  regimes  of  stream  mixing  and  separation. 
Their  experiments  used  photographic  observation  of  the  actual  Impingement  plane 
under  a  transparent  window  and  correlation  of  the  observed  regime  with  the  para¬ 
meters  of  temperature  and  flow  residence  time.  A  review  of  the  problem  was 
presented  by  Kuahkia  and  Houseman  (Ref.  6),  where  an  analytical  model  was 
developed  to  predict  mix/ceparate  criteria  in  terms  of  pressure,  temperature  and 
impingement  plane  contact  time.  This  model  of  liquid  phase  impingement  made 
use  Of  the  heat  release  rate  measurements  of  Feller  and  Somogyi  (Ref.  7),  and 
involved  a  sensible  heat  balance  between  Feller's  heat  release  rate  and  the 
heat  required  to  raise  the  liquid  to  its  boiling  point.  The  model  presented  by 
Law ver  and  Breen  (Ref.  4)  and  that  presented  by  Kushida  and  Houseman  (Ref.  6) 
are  different  in  that  the  first  .'hows  an  exponential  Importance  of  temperature 
upon  the  reaction  r  -te  while  the  second  shows  a  linear  temperature  effect  in 
the  boiling  point  heat  balance,  using  a  maximum  heat  release  rate  from  FeiJer 
and  Somogyi.  It  was  more  the  point  of  view  of  Lawver  and  Breen,  and  is 
expressed  in  Feller's  report,  that  the  hent  release  rate  shows  a  sharp  change 
in  temperature  which  should  correspond  to  tl.e  inception  of  stream  s»  ..aration. 

d.  Summary  of  this  Experimental  Work 

,  The  experimental  techniques  developed  by  Lawver  and  Biaen  in 
Reference  S  were  extended  to  allow  high  pressure  (500  psia)  operation  with  the 
family  of  hydrazine  fuels  N-O.  and  compound  A  oxidizers .  Experimental  data 
wero  thus  obtained  showing  the  effect  of  velocity,  orifice  diameter,  oropellant 
temperature,  impingement  angle  and  chamber  pressure.  These  curve.;  define 
criteria  which  can  be  used  by  the  design  engineer  to  determine  if  hypergolic 
impingement  spray  pattern  is  within  the  mixed  zone  or  the  separate  zone.  At 
longer  residence  times,  popping  phenomenon  was  also  encountered  and  the 
regions  of  popping  are  shown  for  the  hydrazine,  MMH,  and  50-50  fuol^  impinging 
on  N^O..  A  great  many  experimental  difficulties  were  encountered  in  the  obser¬ 
vance  ox  the  compound  A  hydrazine  impingement  plane,  but  It  appeared  that  the 
streams  always  separated  within  the  limited  range  of  temperature  data. 

The  mix/separate  criteria  are  reviewed  in  a  discussion  of  their  applica¬ 
tion.  The  popping  phenomenon,  observed  at  long  residence  time,  has  also  been 
observed  in  engine  test  programs  such  as  that  reported  by  Weiss  and  Gopotek 
(Ref.  3)  and  extensively  investigated  by  Clayton  (Ref.  9).  Beyond  the  immediate 


2 


injector  design,  the  fact  of  Impingement  regime  is  also  discussed  as  it  affects 
such  steady  state  combustion  profile  calculations  as  those  reported  by  Dynamic 
Science  (Ref.  10)  and  currently  being  successfully  used  to  correlate  engine 
stability  by  Abbe  (Ref.  11). 

2 .  IMPINGEMENT  OF  HYPERGOUC  STREAMS 

Both  research  projects  and  full-scale  engine  firings  with  hypergollc 
propellant  combinations  have  demonstrated  that  the  unlike  streams  sometimes 
mix  aud  sometimes  separate  causing  large  changes  in  Injector  and  combustor 
performance.  From  previous  work  it  is  possible  to  develop  a  physical  picture 
of  the  mix/separate  phenomena  and  to  transform  this  picture  or  model  Into 
analytical  expressions  defining  the  controlling  parameters  of  chemical  kinetics, 
residence  time,  and  mixing. 

a  -  The  Physical  Picture 

For  unlike  impinging  hypergollc  streams,  engine  efficiency  depends 
upon  the  degree  of  individual  doublet  mixing  or  separation,  because  generally 
secondary  fan  and  droplet  mixing  is  ineffective-  This  fact  is  demonstrated  by 
the  photographic  comparison  of  the  two  different  injection  schemes  studied  by 
Marshall  Burrows  (Ref.  4).  Figure  la  is  a  6000  frame  per  second  photograph  of 
the  hydrazine  side  of  a  quadruple  impinging  Jet  injector  of  90°  included  angle 
while  Figure  lb  is  the  same  injection  scheme  from  N-O-  side.  The  injection 
scheme  breaks  the  chamber  into  four  quadrants  as  shovm  in  Figures  la  and  lb. 

The  N20.  droplet  vaporization  causes  the  dark  brown  cloud  of  Figure  lb  to  exist 
axially  down  the  chamber;  the  colorless  hydrazine  vapor  of  Figure  la  combusts 
sooner  because  of  the  tendency  for  hydrazine  to  decompose  in  high  temperature 
environments.  In  both  cases  the  combustion  is  efficient  and  is  controlled  by 
droplet  vaporization.  However,  for  the  different  injection  scheme  shown  in 
Figure  lc  the  chamber  is  broken  into  only  two  halfs  and  the  N-O./NLH.  systems 
separate  or  blow  apart.  The  lack  of  O/F  and  mass  distribution  unlfdmrty  causes 
secondary  gas  mixing  to  control  the  combustion  and  the  resulting  stratified  flow 
causes  droplet  vaporization  throughout  the  chamber  to  be  inefficient  because  of 
low  temperatures  and  nonstoichiometric  combustion. 

Observation  of  this  separation  and  stratification  problem  led  to  the 
realization  that  it  would  be  necessary  to  observe  the  details  of  the  phenomena  with¬ 
in  the  impingement  plane  and  to  experimentally  define  the  controlling  parameters. 

A  general  and  complete  description  near  the  stagnation  point  of  impinging 
streams  of  hypergollc  propellants  requires  a  basic  understanding  of  the  inter¬ 
relationships  of  propellant  reaction  kinetics  and  hydrodynamic  flow  properties. 

In  more  detail,  the  reaction  kinetics  involve; 

The  rate  of  liquid  phase  reaction. 

The  rate  of  gas  phase  reaction. 

The  rate  of  heat  release,  and 

The  rate  of  gas  evolution. 

Aside  from  the  stability  problem  associated  with  the  fluid  flow,  important  kine¬ 
matic  parameters  are; 

The  degree  of  mixing. 

The  characteristic  mixing  time,  and 

The  residence  time. 
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End  view  of  inrmgemont  and  resuming 
chamber  quadrants 


Figure  >e.  Efficient  N^O^/N^H^  Combustion  from  the  Hydrazine  Side. 


chamber  qu'drar.ts 

Figure  lb.  The  sens**  Cossbuftor  as  in  Figure  la  showing  Efficient  N204/N2H  . 
Combustion  from  the  Nitrogen  Tetroxlde  Side. 


chamber  quadrants 

Figure  Ic.  Inefficient  N-O./N^H,  Combustion  and  Stratification  Caused  by 
Hypergollo  Stream  Impingement  and  Blow  Apart. 

Fir ura  i.  Tha  Effect  of  Implnqement  on  Stream  Blow  Apart  and  Combustor  Efficiency. 
(Ref  13.) 
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The  condition  of  jet  separation  can  be  viewed  aa  the  occurrence  of  a  ' 
locally  strong  exothermic  reaction/  liberating  sufficient  heat  and  gas  that  the 
streams  are  blown  apart  before  substantial  mixing  has  taken  place.  This 
locally  strong  exothermic  reaction  is  more  likely  to  take  place  close  to  the 
stagnation  point  for  it  is  there  the  propellants  experience  die  longest  residence 
time.  Thus  when  Initial  surface  mixing  first  occurs  near  tha  stagnation  point 
and  the  chemical  reaction  rate  is  fast  enough  to  liberate  sufficient  heat  and 
gas  before  substantial  mixing  has  taken  place  downstream  of  the  stagnation 
point  separated  flow  results.  When  reaction  is  slow  enough  to  allow  a  certain 
degree  of  mixing  to  first  take  place,  a  uniform  combustion  spray  pattern  appears. 
This  physical  picture  has  developed  from  experimental  observation  of  the  impinge¬ 
ment  plane  during  research  performed  on  contract  NA87-467  (Ref.  5).  Typical 
separated  and  mixed  flow  patterns  are  shown  in  Figure  2a  and  2b.  The  propellant 
systems  are  hydrazine  and  nitrogen  tetroxide. 


b.  Data  Correlation 

The  criteria  of  stream  separation  has  been  described  in  detail  in 
Reference  5  as  a  comparison  of  two  characteristic  times,  namely  the  chemical 
kinetic  time  constant,  r  .  ,  and  the  common  interfacial  residence  time  between 

the  two  streams,  r  .  critical  chemical  reaction  time,  r  ,  ,  is  the  time 

required  for  the  chemical  reaction  to  liberate  sufficient  heat  ado^roduct  gas  to 
cause  the  streams  to  separate.  This  time  constant  is  very  sensitive  to  temperature 
because  by  the  very  definition  of  hypergolicity,  a  ’runaway"  ignition  temperature 
exists  within  the  range  of  ambient  propellant  temperatures.  The  Interfacial  resi¬ 
dence  time,  r  ,  is  characterized  by  the  geometry  and  velocity  of  the  flow.  Se¬ 
paration  or  mfimig  phenomena  depend  on  whether  T  .  is  either  smaller  (separate) 
or  greater  (mix)  thanr  .  The  limit  for  chemical  reaction ,  blow-apart,  and 
separation  is  defined  Wen  Treg  Is  equal  to  rc^imm 

The  chemical  kinetic  time  constant  may  be  obtained  by  approximating 
the  rate  of  consumption  of  one  of  the  propellant  species  by  the  first  order  decay: 


dY 

dt 


-kY 


(i) 


where  Y  is  toe  species  concentration  and  k  is  the  kinetic  rate  constant  of  the 
chemical  reaction.  Denoting  Y  as  the  initial  propellant  concentre  lion,  equation 
(1)  gives:  ° 


Y  _  -kt 


(2) 


-1 


and  characteristic  time  to  Ignition,  t  may  be  defined  when  Y/Yq  *  e* 

For  temperatures  below  the  ignition  temperature  this  time  is  very  long,  while  for 
temperatures  greater  than  the  ignition  temperature  the  time  constant  is  much  less 
than  the  characteristic  mass  transfer  on  heat  transfer  time  constants.  Thus  for 
ignition  phenomena  to  occur  the  value  of  the  characteristic  kinetic  time  constant 
shifts  from  rate -controlling  (at  low  temperature)  to  not  controlling  oyer  a  small 
increment  of  temperature  change  and  thus  the  ignition  temperature  Is  closely 
defined  for  any  chosen  value  of  Y/Y  of  the  order  of  one.  Thus  r  .  is  c  ined 
from  equation  (2)  as 


chm 
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The  reaction  rate  constant  k  may  be  assumed  to  follow  the  Arrhenius  temperature 
dependence  so  that 


cltm 


1  ttE/RT 
A  * 


Thus  the  characteristic  time  of  chemical  reaction  is  exponentially  dependent 
upon  temperature. 


(1) 


The  characteristic  residence  time  for  which  the  streams  are  In  forced 
mixing  contact  is  a  function  of  stream  dimensions,  velocity,  and  impingement 
angle.  Characteristically  this  may  be  obtained  from  the  geometry  of  the  situa¬ 
tion  as  ? 

p  i_ 

Tres  ”  sin (6  7 2)  V 

where  0  is  the  Included  angle  of  Impingement,  D  is  the  orifice  diameter,  and 
V  is  the  injected  stream  velocity. 

Whenever  the  chemical  time  constant  is  less  than  or  equal  to  the 
residence  time  the  streams  will  tend  to  blow  apart  due  to  the  tremendous  heat 
release  rates  associated  with  the  ignition  and  combustion  of  stoichiometric 
rocket  propellants  within  the  residence  interface.  Thus  the  separation  limit 
is  defined  in  terms  of  operating  temperature  and  flow  parameters  by  equating 
the  two  time  constants 


E/RT 

e _ 


so  that 

A 

~  sin  (^2)V 

(6) 
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c  sinCe/a)  1  cE/RT 

A 
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or 
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(8) 

Postulation  of  the  impingement  process  in  terms  of  the  physical  parameters 
controlling  reaction  and  stream  separation  indicates  that  the  critical  separation 
limit  varies  with  D/V  and  is  exponentially  dependent  upon  1/T.  Thus  experimental 
data  are  correlated  on  a  semllogarithmic  scale  as  proportional  to  the  log  D/V  and 
Inversely  proportional  to  propellant  temperature. 

c.  Injection  Mixing  Explosions 

Another  very  critical  problem  associated  with  the  impingement  of 
hypergoiic  streams  is  concerned  with  the  often  observed  explosions  near  the 
injector  plate  since  they  can  act  as  initiation  sources  of  detonation  waves  or  ’pops". 
Typical  explosions  of  this  kind  were  Identified  experimentally  at  Dynamic  Science 
under  NASA  research  contract  NAS7-467  (Ref.  12)  and  are  shown  photographically 
In  Figure  3.  In  contrast  to  tne  separation  phenomena  which  are  basically  due  to 
a  locally  strong  exothermic  reaction,  explosions  appear  to  be  caused  by  ignition 
within  the  entire  impingement  region  of  mixed  liquid  ligaments.  In  this  respect. 


Figure  2a.  Mixing  and  Uniform  Spray  Combustion  of  N-O./to-H.  at  14°C. 
at  50  ft/sec  and  .050  Inch  Wide  Stream.  ’  *  * 

Regime  I- Mixing 


Figure  2b.  Hydrazine /Nitrogen  Tetroxlde  Impingement  25°C 
at  lOOft/sec  and  .100  inch  Wide  Stream. 

Regime  II- Separation 


N2H4  VELOCITY  »  25  ft/sec  N^  ORIFICE  DIA  =  .050" 

N204  VELOCITY  »  20  ft/sec  NgO.  ORIFICE  DLA  -  .050" 

PROPELLANT  TEMPERATURE  50  °F 
25  ysec  EXPOSURE  AT  2000  pic/sec 

Figure  3.  Injection  Mixing  Explosions  Causing  Ligament  Shattering  and 
Sometimes  Spray  Detonations  (Taken  from  Reference  12). 


(3)  LIGAMENT  SHATTERING 
INJECTION  MIXING  EXPLOSION 
CAUSING  LIGAMENT  SHATTERING 


NOT  REPRODUCIBLE 


(2)  INJECTION  MIXING  EXPLOSION 


(31  SPRAY  DETONATION 
INJECTION  MIXING  EXFLOSION 
CAUSING  SPRAY  DETONATION 


(2)  INJECTION  MIXING  EXPLOSION 
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Injection  mixing  explosions  depend  very  much  on  the  degree  of  mixing  of  the 
Impingement  region.  First  of  all  the  temperature  must  be  low  enough  that 
mixing  phenomena  predominate  then  the  degree  of  mixing  of  Impinging  Jets 
depends  on  contact  surface  area,  intrinsic  turbulent  levels  and  turbulent 
shear  stresses.  Extensive  experimental  work  (Ref.  I)  has  been  carried  out 
to  measure  the  mixing  level  of  nonreactive  impinging  streams.  Maximum 
mixing  is  observed  when  the  quantity  PiUi  Di  *  paUlDa  -  P  and  U  are  the 
stream  density  and  velocity.  D  is  the  orifice  diameter.  Based  on  experimental 
findings  degree  of  mixing  can  be  defined  in  terms  of  a  ndndimenslonai  number  o  . 


fi  » 


1  + 


Pi1*!0! 

P^ 


(9) 


Uniform  mixing  occurs  at  6  =  0.5.  It  is  anticipated  that  for  6  close  ‘.oO.S, 
injection  mixing  explosions  are  more  likely  to  occur.  Figure  4  shows  ragions 
of  stream  mixing  and  separation  as  well  as  injection  mixing  explosions  as 
functions  of  propellant  temperature,  D/V,  and  degree  of  mixing  arbitrarily 
defined  in  terms  of  6. 

The  present  experiment  is  designed  to  study  stream  n  Ix/separate 
phenomena.  Propellant  temperature  T  and  the  mixing  time  D/V  were  accurately 
varied  in  order  to  Investigate  the  relationship  exoressed  in  equation  (8). 
Photographic  observations  of  the  impingement  region  were  recorded  to  define 
the  mix/separate  phenomena.  The  measurement  of  injection  mixing  explosions 
was  not  the  orimary  concern  of  this  research  but  could  be  observed  qualitatively 
from  chamber  pressure  noise  and  roughness. 
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3.  EXPERIMENTAL  APPARATUS 

®.  Experimental  Technique 

The  experimental  technique  lc  bated  rn  photographic  observation 
of  the  Impingement  region  of  a  pair  of  rectangular  jets.  Injector  plates  of 
different  orifice  sizes  and  impingement  angles  are  used.  Various  Injector 
plates  are  shown  In  Figure  5.  A  pyrex  or  calcium  fluoride  window  which 
extends  beyond  the  Impingement  point,  covers  the  front  face  of  die  Injector 
plate.  This  allows  photographic  observations  of  the  Impinging  streams 
flowing  over  the  window  and  at  the  same  time  still  retains  the  basic  feature 
cf  a  set  of  three  dimensional  impinging  jets.  Liquid  interface  of  the  Impinging 
fuel  and  oxidizer  streams  can  be  distinctly  observed.  technique  allows 

liquid  Interface  of  the  Impinging  streams  to  be  clearly  observed  as  compared 
with  the  Impingement  of  two  free  streams.  Mixing  and  separation  phenomena  are 
identified  and  defined  based  on  these  photographic  results. 

b.  Impingement  Model 

The  stainless  steel  Impingement  test  model  Is  shown  in  Figure  6* 

The  Injector  support  and  the  orifice  plate  are  fitted  into  a  chamber  block. 

The  chamber  is  open  along  Its  whole  length  on  one  side.  A  luclte  window 
with  pyrex  or  calcium  fluoride  Inserts  is  mounted  on  the  open  side  of  the 
chamber.  Ihls  allows  photographic  observation  of  the  Impingement  region. 

A  nozzle  Is  attached  to  the  bottom  of  the  chamber  block.  Chamber  pressures 
of  500  psi  can  be  easily  attained  by  varying  the  nozzle  sizes.  A  schematic 
of  the  chamber  block  and  the  Injector  support  is  shown  In  Figure  7.  The 
support  is  designed  for  a  triplet  as  well  as  a  doublet  Injector  plate.  The 
propellants  are  temperature  controlled  up  to  the  injector  plate  by  means  of 
coaxial  lines. 

c.  Flow  System 

A  schematic  of  the  nitrogen  tetroxlde/hydrazlne  flow  system  Is 
shown  In  Figure  8.  The  flow  systems  are  fabricated  with  304  stainless 
steel  lines.  The  propellant  tanks  are  designed  for  a  working  pressure  of 
1000  psi.  The  tank  capacities  are  approximately  1/3  gallon.  Pressurization 
of  the  tanka  is  controlled  by  a  Grove  regulator  and  a  release  valve  which 
Is  set  at  1000  psi.  Two  solenoid  valves  are  used  to  vary  the  dome  pressure 
of  the  regulator.  The  dome  press;.’’’®  is  always  set  at  approximately  950  psi. 
Propellant  flow  rates  are  controlled  by  variable  area  cevitating  venturiea  and 
are  metered  by  turbine  type  flow  meters.  The  flow  controller*  are  micrometer 
head  type.  This  device  provides  easy  and  precise  control  of  the  propellant 
flows  and  allows  for  more  efficient  testing  with  the  ability  to  match  desired 
stream  momentum  ratios.  The  propellant  tanks  and  lines  are  provided  with 
heat  exchangers.  Propellant  temperature  is  controlled  by  flowing  coolant 
at  the  desired  temperature  through  the  heat  exchangers.  The  coolant 
temperature  Is  regulated  with  a  refrigeration  unit  with  heat  exchangers  on 
the  condenser  and  evaporator.  Coolant  temperature  is  controlled  by  throttling 
the  flews  through  the  condenser  and  evaporator  heat  exchangers.  The  propel¬ 
lants  are  first  led  lo  flow  to  the  exhaust  tank.  When  the  desired  temperature 
is  achieved,  xho  propellant  is  introduced  into  the  impingement  block.  Pres¬ 
surization  of  the  exhaust  tank  is  achieved  by  a  two  stage  regulator.  This 
enables  recovery  of  the  propellants  from  the  exhaust  tank. 
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Flow*  7a .  Injector  Mounting  Support  for  Hot  Firing 


Figur*  7b.  Impinging  Elocfc  for  Hot  Firing 


A  separate  flow  system  is  used  to  provide  liquid  compound  A  to  the 
impingement  apparatus.  The  schematic  of  the  liquid  compound  A  supply 
system  is  shown  in  Figure  9.  The  system  is  designed  for  a  working 
pressure  of  1000  pal.  Compound  A  temperature  can  be  easily  varied  from 
-300°F  to  120°F.  All  the  nitrogen  refrigerant  and  purging  gas  is  provided 
by  an  809  gallon  liquid  nitrogen  storage  dewar  maintained  at  40  psl,  snd 
located  adjacent  to  the  test  area.  This  supports  the  test  setup  with  extremely 
dry  nitrogen  gas  for  the  various  purging  cycles  and  cold  liquid  and  gas 
nitrogen  to  provide  low  temperature  conditioning.  Propellant  temperatures 
above  ambient  conditions  is  achieved  by  heating  the  nitrogen  gas.  Compound 
A  is  stored  in  cylinders  and  transferred  to  the  run  tank,  as  required,  during 
the  test  program.  Pressurization  of  the  run  tank  is  achieved  with  helium  pcs. 
To  insure  tht  moisture  content  of  the  helium  gas  is  within  acceptable  limits, 
ali  helium  used  for  propellant  press u-ization  is  passed  through  a  liquid 
nitrogen  refrigerated  molecular  sieve  dryer. 

The  vent  gas  of  the  flow  system  is  neutralized  with  a  charcoal 
reactor  before  being  exhausted  into  the  atmosphere.  Combustion  products 
from  the  lest  chamber  are  aspirated  with  nitrogen  gas,  scrubbed  with  water, 
and  periodically  neutralized  in  the  disposal  ponds. 

The  flow  system  is  being  operated  remotely  during  tent .ins.  Various 
solenoid  valves  and  pneumatic  controlled  ball  valves  are  used  for  the  hydra¬ 
zine  and  nitrogen  tetroxide  supply  systems.  Pneumatic  controlled  Ni/pro 
bellow  valves  are  used  for  the  compound  A  system. 

Pressure  is  measured  by  strain  gauge  transducers.  The  strain  gauge 
is  powered  by  a  six-volt  power  supply.  The  transducer  output  is  fed  through 
the  patch  panel  into  an  amplifier  to  the  oscillograph  galvanometer.  The 
strain  gauge  system  Is  provided  with  a  zero  balance  circuit  and  an  electrical 
calibration  unit  which  is  put  l’-to  the  circuit  through  a  selector  switch.  The 
pressure  transducers  have  a  natural  frequency  of  8000  CPS  and  a  linearity 
of  +  0.5%.  The  transducers  are  calibrated  with  nitrogen  gas  using  a  “dgh 
accuracy  Bourdan  gage.  The  Bourdan  gage  is  calibrated  to  an  NBS  standard 
and  is  accurate  to  +  0 . 25%. 

Hydrazine  and  nitrogen  tetroxide  flow  rates  are  measure J  wi;.h  Cox 
turbine  type  flow  meters.  Compound  A  flow  rates  are  measured  with  two 
Potter  flow  meters.  The  flow  meter  output  signals  are  fed  into  an  amplifier 
to  the  oscillograph.  The  flow  rate  is  determined  from  the  meter  frequency 
output  and  the  propellant  density. 

Propellant  temperatures  are  measured  with  I/C  thermocouples.  The 
thermocouple  output  is  fed  through  the  patch  board  into  an  amplifier  to 

the  oscillograph. 

The  impingement  test  stand  with  compound  A  supply  system  is  3hown  In 
Figure  IQg.  A  close  view  of  the  test  stand  is  shown  in  Figure  10b, 
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Figure  9.  ClFj.  Supply  System 
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d.  Control  Console 

Tho  test  control  console  houses  the  facility  and  propellant  control 
panel,  a  test  sequence  unit,  amplifier,  and  recording  oscillograph  with  related 
components.  The  console  is  shown  in  Figure  1?  The  Instrumentation  system 
Is  centered  around  a  Midwestern  801  recording  oscillograph.  Propellant 
temperatures  at  the  injector  plates,  propellant  flow  rates  and  the  chamber 
pressure  are  recorded  on  the  oscillograph.  Propellant  tank  pressures  are 
recorded  on  a  two  channel  strip  recorder.  Propellent  temperatures  downstream 
of  cavltatlng  venturies,  tank  temperatures,  and  flow  meter  temperatures  are 
recorded  on  a  multipoint  recorder. 

e.  Photography 

The  Impingement  region  Is  photographed  using  a  Graphex  4x5  still 
camera.  A  Written  No.  34  filter  is  used  to  mask  out  the  flame  light.  This 
filter  transmits  only  in  the  blue  and  red  region.  Polorold  color  film  Is  used 
with  a  carbon  arc  lamp  and  photoflood  front  lighting  techniques . 
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Figure  11.  Control  Console  with  Facility  Control  Panel. 
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4.  IMPINGEMENT  OF  NITROGEN  TETROXIDE/ HYDPAZIN^  FAMILY 

PROPRLLANT  SYSTEMS. 

a.  Experimental  Data 

Impingement  experiments  were  conduc  ed  by  Impinging  streams  of 
hydrazine,  monomethyl -hydrazine  and  mixtures  of  hydrazine  and  unsymmetric- 
di-methyl  hydrazine  (A-50)  with  nitrogen  tetroxld*.  Propellant  temperatures  were 
varied  from  40°F  to  120°F.  The  characteristic  mixing  time,  D/V  was  varied 

-4  -4 

from  0.2x10  seconds  to  3x10  seconds.  To  evaluate  the  effect  of  pressure,  * 
chamber  pressure  is  varied  between  15  psi  to  500  psf. 

Figures  12,  13,  and  14  are  photographs  of  impingement  region  of  hydrazine 
family  fuels  with  nitrogen  tetroxide.  A  Wri'.ten  No.  34  filter  is  used  to  mask  out 
the  flame  light.  A  dark  purple  region  along  the  contact  surface  is  present  for 
monomethyl-hydrazlne/nltrogen  tetroxide  impingement  (Fig.  13a).  This  is  most 
probably  due  to  certain  hydrocarbon  intermediates. 

Figure  15  is  the  impingement  data  for  nitrogen  tetroxide  and  hydrazine 
streams  (Ref.  5)  for  both  rectangular  ai.d  circular  Jets.  The  data  are  presented 
in  terms  of  the  propellant  temperature  expressed  as  1/T  and  the  mixing  time 
D/V.  The  experiments  are  conducted  at  the  ambient  pressure  of  15  psi.  The 
experimental  data  in  terms  of  the  Inverse  of  the  absolute  propellant  temperature 
and  D/V  agree  very  well  with  the  relationship  suggested  in  Equation  (8).  Figure 
16  is  nitrogen  tetroxide/hydrazL  e  impingement  data  for  chamber  pressures  varying 
between  15  psi  and  500  psi.  The  angle  of  Impingement  Is  generally  60.  Impinge¬ 
ment  angles  of  30°  and  £0°  are  also  used  occasionally.  Chamber  pressure  does 
not  seem  to  play  any  significant  role  with  respect  to  stream  separation  phenomena. 
Fioures  17  and  18  are  Impingement  data  for  nitrogen  tetroxlde/monomethyl  hydra¬ 
zine  and  nitrogen  tetroxide/aerozlne  50  respectively.  Separation  limit  for  these 
propellants  differs  significantly  from  that  of  nitrogen  tetroxide  and  hydrazine. 

b.  Discussion 


The  impingement  model  outlined  in  Section  2b  relating  the  propellant 
temperature  and  the  mixing  time  (Eqn.  8)  allows  good  correlation  with  experimental 
data.  Beyond  this  correlation  the  order  of  magnitude  of  the  "overall  activation  energy" 
E,  car.  be  estimated  based  on  the  slope  of  the  limiting  curve.  F.  assumes  a  value  of 
approximately  20  keal/moie  of  reactants.  A  more  detailed  calculation  has  been 
carried  out  by  Jaye  and  Blauer  (Ref.  14)  for  the  date  of  Figure  15.  By  assuming 
that  Ignition  occurs  at  1%  of  the  total  reaction  and  using  sin  6/2»l/2  for  a  6(ftmpinge- 
ment,  >ey  calculated  that  E  and  log  A  are  respectively  24  +  3  kcal/mole  and 
15.  /  0.5  cc/mole/sec. 

As  observed  by  Jaye  and  Blat.r,  rate  constants  obtained  from  the  Impinge¬ 
ment  results  are  m  excellent  agreement  with  the  work  of  Sawyer  and  Glassman, 

{Ref.  15).  They  obtained  E  -  26.7  +  1.1  kcal/mole  and  log  A«15.8  +  0.o  cc/mole/sec 
for  the  gas  phase  Ignition  process.  Jaye  and  BJauer  continue  that  this  agreement 
leads  "to  the  conclusion  that  either  the  rate  controlling  process  (of  stream 
separation)  occurs  in  the  gas  phase  or  the  reaction  mechanism  is  the  same  in  both 
the  liquid  and  gas  phases."  This  very  interesting  observation  is  supported  by 
consideration  of  surface  physics.  When  two  liquid  surfaces  are  In  contact,  the 
surface  layer  molecules  which  contact  the  opposite  surface  would  tend  to  follow 
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Figure  12a.  Stream  Doublet  Separation  (N2H/,/N204). 


Figure  12b.  Stream  Doublet  Mixing  (N~H./N_0,h 
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Figure  13*  Stream  Doublet  Separation  (MMH/N^). 


Figure  13b.  Stream  Doublet  Mixing  '(MMH/N.,0^). 
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Figure  i  6 .  Separation  Data  for  ^O./N^H.  System  for  Different 
Chamber  Pressures. 
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vapor  pressure  kinetics  rather  than  fully  mixed  liquid  phase  kinetics.  This 
model  is  supported  by  motion  pictures  of  Weiss,  Fisher,  Gersteln,  and  Johnson 
{r«.er.  i  6)  through  which  they  observed  that  nitrogen  tetroxlde  drops  appeared 
to  be  immissible  in  hydrazine  and  maintained  their  original  surface  Identify. 

The  data  presented  in  Figure  16  definitely  indicate  that  stream  separation  does 
not  depend  on  the  chamber  pressure.  When  a  locally  strong  exothermic 
reaction  takes  place  between  die  propellants  to  liberate  sufficient  gas  pro¬ 
ducts  and  heat  to  vaporize  the  propellants  and  the  reaction  intermediates  so 
as  to  blow  the  stream  apart,  chamber  pressure  only  plays  a  secondary  role  of 
increasing  the  boiling  point  of  the  liquid.  Since  the  latent  heat  is  the  dominat¬ 
ing  factor,  the  effect  of  chamber  pressure  does  not  play  a  significant  role.  To 
further  illustrate  this  point,  the  boiling  point  of  NjO,  at  IS  psi  and  500  psl  are  530  R 
and  700°R  respectively.  The  latent  heat  of  vaporization  within  these  temperature 
ranges  are  approximately  400  Btu/lb  and  the  specific  heat  C  Is  approximately 
0.3  Btu/lb°R.  The  significance  of  pressure  can  be  evaluated  by  the  ratio: 


Cn[Th(500  psl)  -  Tb(15  psl)]  0.3(1701  „„ 

latent  heat  "  400  “ 

This  indicates  that  a  S00  psl  pressure  increase  only  causes  a  first  order  effect  on 
the  separate/mix  temperature  criteria.  When  the  order  of  magnitude  cf  the  ratio  of 
the  sensible  heat  difference  to  the  combustion  heat  available  is-  compared,  the 
amount  of  heat  used  in  reaching  different  boiling  points  is  even  less  significant. 

Preliminary  tests  on  nitrogen  tetraxide/aerozine  50  and  nitrogen  tetroxlde/ 
monomethyl-hydrazlne  systems  have  been  carried  out.  Separation  limits  of  these 
propellant  systems  are  to  be  compared  with  that  of  nitrogen  tetroxide/hydrazine 
systems.  For  a  given  mixing  time  D/V,  separation  limits  of  these  propellant 
systems  occur  at  a  higher  temperature  than  that  of  hydrazine.  This  indicates  that 
the  nitrogen  tetroxide/hydrazine  is  a  more  reactive  propellant  system  as  compared 
with  that  of  nitrogen  tetroxlde/monomethyl-hydrazlne  and  nitrogen  tetroxlde/ 
aerozine  50  systems.  Based  on  jet  separation  phenomena,  at  a  given  temperature 
and  mixing  time,  aerozine  50  and  monomethyl -hydrazine  with  tetroxlde  will  more 
readily  mix  and  therefore  give  a  more  intense  initial  combustion  zone  than  the 
hydrazine  with  nitrogen  tetroxlde.  However  In  each  case  the  separation  temperature 
may  be  reached,  depending  upon  the  use  of  regenerative  cooling  and  resulting 
injection  temperature. 

c.  Engine  Design  Application 

Stream  separation  influences  the  engine  performance  In  many  respects. 
When  reactive  streams  are  blown  apart  the  result  can  be  a  high  degree  of  flow 
striation.  This  leads  tc  poor  combustion  efficiency.  When  separation  occurs 
maximum  striation  and  efficiency  loss  is  associated  with  a  single  element  injector. 
The  striation  effect  upon  efficiency  losses  would  diminish  for  larger  numbers  of 
injector  elements.  This  effect  was  established  experimentally,  (Ref.  17),  and 
indicated  that  the  combustion  efficiency  loss  for  a  single  element,  four  on  one, 
injector  Is  roughly  twice  that  of  a  ten  element  injector. 

The  influence  of  nonuniform  mixture  ratio  distribution,  which  can  be 
induced  by  stream  blow  apart,  has  been  theoretically  calculated  based  on  a 
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stra titled  supersonic  nozzie  flow  field  (Ret .  18)  and  internally,  based  on  non- 
uniform  vaporization  streamtubes  (Ref.  19).  Varioue  prescribed  oxidizer  to 
fuel  ratio  profiles  havo  bean  used  to  evaluate  their  affect  on  performance. 

Nitrogen  tetroxlde/aerozine  SO  and  fluoriny/hydrogen  systems  generally  indicated 
; arses  from  1.0  to  4.0%  in  specific  impulse.  Liquid  oxygen/h^drocarbon 
(LCX'RP-J)  system  exhibited  a  specific  Impulae  loss  of  approximately  8%. 

The  mix/ separate  criteria  presented  in  this  report  also  point  out  the 
importance  of  proper  injector  design  upon  thrust  chamber  liner  compatibility. 

Proper  injector  design  must  include  the  Influence  of  propellant  temperature 
upon  whether  or  not  the  fuel -rich  stream  separates  and  remains  next  to  the  liner; 
or  in  the  low  temperature  case  (mixing)  forms  an  intfeuse  combustion  zone.  This 
effect  will  be  found  not  only  in  the  fuel -rich  barrier  zone  but  also  in  the  Intense 
stoichiometric  combustion  of  the  primary  zone.  If  the  primary  injector  zone  is 
within  the  mixing  fiegime  it  may  actually  blow  through  the  fuel-rich  combustion 
of  the  barrier  -one.  Thus  liner  compatibility  is  as  strongly  dependent  upon 
propellant  temperature  as  it  is  upon  propellant  dist 'button. 

It  has  bean  shown  in  the  previous  section  that  by  properly  controlling  the 
propellent  temperature,  propellant  injection  velocity  and  orifice  size,  regimes 
of  mlx/separate  phenomena  o«m  be  defined.  This  knowledge  can  be  used  to 
great  advantage  by  design  engineers.  Based  on  data  similar  to  those  presented 
in  Figures  IS  to  18,  characteristic  parameters  such  as  propellant  temperature 
and  mixing  time,  D/V  can  be  defined  in  order  to  achieve  proper  combustion 
mixing  leading  to  better  engir.~  performance,  stability,  and  compatibility.  Parti¬ 
cularly  when  Anomalous  engine  behavior  is  observed,  attention  should  be  paid 
to  the  oper  t4ng  region  of  mix  or  separate  mechanisms,  et  the  very  least  this  work 
points  out  the  Importance  of  controlling  and  monitoring  propellant  temperature  during 
engine  firing  test  programs. 

Engine  instability  is  another  problem  which  is  sometimes  directly  attributed 
to  separation  phenomena.  High  pressure  transients  have  often  been  observed  at 
chamber  wall  pressure  transducer  stations  during  full-scale  engine  development  and 
test  programs  (Refs.  3  and  9).  These  transients  are  rather  arbitrarily  termed  "pops," 
depending  upon  amplitude  and  frequency.  However,  their  initiation  source  is  the 
subject  of  much  speculation.  During  the  course  of  this  Impingement  research  and 
in  previous  impingement  experiments  conducted  at  NASA/Lewis  (Ref.  13)  and  at 
Dynamic  Science  (Ref.  12),  it  was  found  that  injection  mixing  explosions  often 
occurred  teproduclbly  in  die  lower  temperature  hr  lx)  regime,  at  long  residence  time. 
Their  occurrence  also  appears  co  increase  under  these  conditions  when  the  streams 
besvs  a  high  degree  of  mixing  as  shown  in  Figure  4,  (i , e .  5  defined  on  page  7  to  be 
approximately  0.5  or  that  the  streams  have  approximately  equal  momentum  or 
stagnation  pressures,) 

High  speed  motion  picture  cameras  have  been  used  to  record  the  explosion 
of  the  impingement  region  (Ref.  .12).  This  explosion  phenomena  is  very  detrimental 
to  stable  engine  operation  because  they  can  serve  as  triggering  source  for  spray 
detonation  as  shown  in  Figure  3,  This  detonation  wave  at  the  impingement  zone 
also  causes  Injector  flow  irregularities  with  apparent  "flow  shutoff"  of  u  few  milli¬ 
seconds  resulting  in  rough  engine  operation.  The  pictures  on  Figure  3  were  obtained 
with  long  interface?  residence  time  and  substantial  liquid  mixing  throughout  the 
impingement  region.  Unlike  the  separation  phenomena  which  are  caused  by  a 
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locally  exothermic  reaction,  the  reaction  within  the  entire  impingement  region 
resultc  in  spray  detonation.  For  smooth  engine  operation  the  design  engineer 
has  to  avoid  these  conditions  which  lead  to  detonation  waves.  The  occurrence 
of  these  explosions  will  be  unavoidable  with  unlike  hypergolic  propellant  injec¬ 
tion  unless  care  is  taken  to  insure  that  each  set  of  elements  within  the  injector 
are  operating  within  the  same  regime.  Because  of  manifold  layout  and  I/D 
orifice  problems  this  is  often  impossible  but  it  should  always  be  a  major 
consideration  for  the  removal  of  "pop'"  and  roughness.  If  each  set  of  elements 
are  different  there  may  be  a  random  sc:  sometimes  passing  through  the  Injector 
mixing  explosion  regime.  The  use  of  a  single  element  system  such  as  a  center 
flow  injection  scheme  overcomes  the  element  randomness. 
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S .  COMPOUND  A  IMPINGEMENT  WITH  HYDRAZINE 
a.  Experimental  Difficulty* 

The  Impingement  model  (Fig.  G)  for  nitrogen  tetroxide/hydrazlne 
taste  Is  also  used  to  stud/  the  jet  separation  of  compound  A  and  hydrazine. 

A  separate  flow  system  (Figs.  7  and  8a)  is  used  for  delivering  compound  A. 
Propellant  is  temperature  regulated  up  to  the  injector  plate  (Figure  7a), 
since  heat  transfer  between  the  plate  which  is  initially  at  ambient  tempera¬ 
tures  and  compound  A  which  is  close  of  0°F  only  increases  the  propellant 
temperature  by  a  fraction  of  a  degree. 

Impinging  compound  A  with  hydrazine  presents  many  additional  diffi¬ 
culties.  Compound  A  is  a  more  reactive  and  corrosive  propellant  as  compared 
with  nitrogen  tetromide.  This  often  causes  substantial  damage  to  the  viewing 
window  of  the  Impingement  test  model.  Luclte  is  certainly  not  compatible  with 
compound  A  and  this  is  evidenced  with  die  initial  runs  using  lucite  windows. 
Reactions  between  luclte  and  compound  A  upstream  of  the  impingement  point 
is  shown  in  Figure  19.  .  The  bright  flame  along  the  entire  length  of  the  compound 
A  nozzle  is  the  reaction  of  compound  A  with  lucite.  This  in  turn  caused  sub¬ 
stantial  damage  to  the  injector  plate  as  shown  in  Figure  20.  To  partly  remedy 
this  problem,  the  injector  plate  lr  modified  as  shown  in  Figure  21.  This  also 
allows  impingement  to  take  place  away  from  die  window.  Luclte  with  pyrex 
or  calcium  fluoride  inserts  are  used  for  window  materials.  Damages  on  these 
windows  are  very  often  too  severe  to  obtain  quality  photographs.  The  chamber 
window  following  a  compound  A/hydrazine  impingement  reaction  is  shown  in 
Figure  22.  It  is  a  lucite  window  with  a  pyrex  Insert  1/4  inch  thick.  In  addi¬ 
tion  to  the  severe  bumhg  of  the  window  by  compound  A,  much  damage  is  caused 
by  the  thermal  stresses  resulting  from  the  strong  exothermic  reactions  at  the 
Impingement  region.  A  typical  photograph  using  calcium  fluoride  Insert  is 
shown  in  Figure  23.  By  closely  investigating  the  calcium  fluoride  insert  following 
t*e  Impingement  test,  damage  was  found  to  be  mainly  caused  by  thermal  stresses 
since  the  window  is  3/8  inch  thick.  It  appears  that  calcium  fluoride  is  com¬ 
patible  with  the  poopound  Aat  these  high  temperatures.  Better  quality  photo¬ 
graphs  might  be  obtained  with  thin  calcium  fluoride  inserts. 

The  vent  gas  of  the  flow  system  is  neutralized  with  a  charcoal  reactor 
before  vented  into  the  atmosphere.  Refractory  cement  is  used  as  e  heat 
insulator  for  the  charcoal  scrubber.  The  strong  exothermic  reaction  between 
compound  A  and  charcoal  sometimes  'melts*4  this  cement  Insulator  which  is 
shown  in  Figure  24. 

b.  Results 

A  photograph  of  impinging  streams  of  hydrazine  and  compound  A  is 
shown  in  Figure  25.  The  fuel  and  oxidizer  flow  rates  are  0.022  Ib/sec  and 
0.0343  iiV'sec.  The  initial  temperatures  are  56°F  and  0SF.  based  on  photographs 
of  the  impingement  region,  mLx/separate  phenomena  is  less  distinct  as  compared 
with  nitrogen  tctrcxxide  and  hydrazine.  Unlike  nitrogen  tetroxide  which  is  brown 
colored  Compound  A  is  colorless. 

Approximately  fifty  impingement  tests  using  hydrazine  and  compound  A  have 
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Figure  23.  Luclte  Window  with  Calc'um  Fluoride  Insert. 


Figure  24. 


Alelteu  Refractory  Cement  of  Charcoal  Scrubber. 
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Figure  25.  Hydrazine/Compound  A  Impingement. 
(Compound  A,  Left  Injector 
Hydrazine,  Right  Injector). 
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been  carried  out.  The  chamber  pressure  was  varied  from  ISpsia  to  over  iOO  psia.  m 

Hydrazine  and  Compound  A  temperatures  ranged  from  40°F  to  6D°F  and  0CF  to^  M 

30°F  respectively.  The  mixing  time,  D/V  was  w*thin  0.2xl0~*  sec  to  2x10  sec.  J 

Within  these  temperature  and  velocity  ranges,  it  appeared  that  impingement  of  ji 

compound  A  and  hydrazine  always  resulted  in  separation.  ^ 

The  separation  phenomena  are  partly  attributed  to  the  high  vapo  pressure 
of  compound  A.  Its  boiling  temperature  at  atmospheric  pressure  is  7°F,  wnlle 
the  hydrazine  had  to  be  maintained  above  35°F  to  keep  it  from  freezing.  Thus  « 

the  compound  A  would  boil  on  contact  with  the  hydrazine  at  pressures  below 
45  psia  (corresponding  to  hydrazine  temperature  of  60°F).  Under  such  circum¬ 
stances  chamber  pressure  is  indeed  Important  to,  first  of  all,  assure  liquid- 
liquid  contact.  The  significance  of  chamber  pressure  in  this  case  is  not  con¬ 
tradictory  to  the  data  of  Figure  16  which  show  no  pressure  effect,  but  merely 
indicate  that  a  boiling  mechanism  will  cause  separation  up  to  the  point  where 
the  chamber  pressure  allows  liquid-liquid  contact. 

It  is  believed  that  better  experimental  results  can  be  achieved  with  the 
viewing  window  further  away  from  the  impingement  region.  A  lucite  window  with 
calcium  fluoride  Inserts  was  compatible  with  compound  A  at  *he  elevated  tempera- 
ures.  In  the  present  experiment,  test  periods  are  approximately  two  seconds. 

Reducing  testing  periods  will  minimize  possible  window  damage.  Ihis  can  be 
achieved  by  improving  the  existing  by-pass  line  in  the  flow  system  so  as  to 
aduce  the  transient  effects.  Possibility  of  coloring  one  of  the  propellants 
can  be  investigated  so  as  to  obtain  better  photographs. 
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6.  CONCLUSION 


BaMd  on  experimental  results  of  hypergolir  propellant  impingement 
studies,  mixing  or  separation  phenomena  are  usually  observed .  Corresponding 

A#  first  Damkohlar  number,  a  nood  linen  signalized  number  3  -  r  .  /r 
.  .  .  chm  mix 

can  be  defined.  fchfl|  and  rmU  are  characteristic  chemical  reaction  time  and 

mixing  fan.  Separation  limit  is  defined  by  assuming  3*1. 

Impinging  nitrogen  tetroxlde  with  hydrasine  family  fuels,  the  separation 
limit  3-  1  agrees  very  well  with  the  experiment.  Chamber  pressures  * 
varied  boa  15  pel  to  500  psl.  The  predicted  pressure  effects  (Ref.  6#  ware 
not  observed.  In  fact,  in  the  regime  that  was  investigated,  stream  separation 
iff  independent  of  pressure.  This  is  anticipated  when  chemical  kinetics  are 
go  versed  by  the  rate  of  liquid  phase  reaction.  Separation  limits  are  observed 
at  higher  temperatures  for  A-5Q/N204  and  MMH/N2C4  as  compare.’  '.th 
N2VN2G4  system  indicating  l$2H4/fc204  is  a  mc*a  reactive  system. 

It  appears  shat  stream  separation  always  occurs  by  impinging  Compound  A 
and  hydrazine  with  oxidizer  and  fuel  temperatures  varying  between  0°  F  to  30°F 
and  40*P  to  60*F  respectively.  Better  experimental  results  can  be  obtained 
by  reducing  the  testing  period  and  using  lucite  windows  with  thin  calcium 
fluoride  inserts. 
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PART  II 


EXTENSION  OP  STEADY- STATE  COMBUSTION  MODEL 


i .  IN  TROD  UCTXOM 

The  Dynamic  Science  steady-state  combustion  computer  program  has 
been  in  use  at  several  installations  for  a  number  of  years.  It  has  been  shown 
by  many  Investigators  to  be  a  useful  tool  to  investigate  the  performance  and 
stability  of  liquid  rocket  engines.  The  program  provides  the  linkage  between 
propellant,  engine,  and  injector  parameters,  and  die  performance  of  the  system 
and  guides  the  designer/analyst  in  making  changes  necessary  to  correct  sys¬ 
tem  deficiencies 

This  section  describes  a  project  performed  jointly  by  personnel  of  the 
jJrltPL  and  Dynamic  Science  under  the  direction  of  Capt.  Charles  J.  Abbe  to 
improve  the  operational  versatility  and  economics  of  use  of  the  steady-state 

combustion  program. 


2 .  PROGRAM  DESCRIPTION 

The  basic  model  used  in  the  Steady-Ctate  Combustion  Program  is  an 
extension  of  the  original  work  of  Priem  (Ref.  20).  The  model  considers  the 
conservation  equations  governing  the  heat,  mass,  and  momentum  transfer  of  a 
spray  of  liquid  fuel  drops  traveling  through  a  channel  filled  with  oxidizer  vapor. 
The  assumption  is  made  that  propellant  vaporization  is  the  rate  governing 
mechanism  of  the  combustion  process  and  that  vaporization  of  the  drop  can  be 
expressed  by  a  Ranz  and  Marshall  type  correlation. 

Dynamic  Science  extended  this  work  (Refs.  21  and  22)  to  Include  the 
simultaneous  vaporization  of  fuel  and  oxidizer  sprays;  thereby  calculating 
vapor  and  liquid  phase  O/F  along  the  chamber  length.  The  thermophy3ical 
property  calculations  were  also  substantially  improved,  A  valuable  routine 
add^d  to  the  program  included  the  monopropellant  burning  effect  of  hydrazine 
type  propellants  and  allowed  calculated  and  measured  combustion  efficiencies 
to  be  correlated  with  these  propellant  systems  (Ref.  23). 

The  program  as  listed  in  Reference  20  included  all  these  effects  but 
was  still  in  an  engineering  development  stage  not  entirely  suited  to  production 
running.  This  work,  in  this  computer  extension  phase.  Involved  breaking  the 
program  up  into  subroutines  of  size  suitable  for  manipulation,  replacing  the 
second  order  Runge-Kutta  integration  with  a  predictor-corrector  type  scheme, 
and  improvement  of  program  versatility.  The  bulk  of  the  work  was  performed 
for  the  Improvement  of  program  versatility,  i.e. ,  generalizing  the  program  to 
provide  the  capability  of  easily  switching  from  one  propellant  combination  to 
another.  This  Involved  segregating  the  property  routines  and  generalizing  their 
input,  gathering  the  property  data  available  and  converting  to  the  proper  units 
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sac!  specifying  estimation  procedures  for  properties  not  readily  available.  The 
fisssjorlty  oi  proaramming  and  ail  debugging  were  performed  at  the  AFRt*.  computing 
fac''*ty. 

Thfe  documentation  citod  here  la  only  an  attempt  to  provide  a  knowledge 
ol  the  changes  that  have  been  mad*  tn  the  program .  For  a  more  complete 
description  of  the  modal  and  program,  see  References  13  *nd  20. 

3 .  PHYSICAL  PROPERTY  RO  UTINES 

The  following  Is  a  list  of  the  pnysicai  property  routines  and  uni *5 
required  by  the  program.  V*  *  properties  may  be  broken  down  into  those  of  the 
fuel  and  o..*dlzsr  (both  liquid  and  vapor!  and  those  cf  the  individual  species 
comprising  the  combustion  gss,  ‘£hs  independent  variable  against  whith  values 
are  to  be  Interpolated  is  also  given.  Oxidizer  properties  may  be  relaced  in  the 
program  tc  the  integer  I  *  1  end  fuel  properties  by  I  *  2. 


_ _  _ _ _  Fuel  and  Oxidizer  Properties _ _ 

Subroutine  Property  (Units)  Independent  Variable  (Units) 

FVAP  vapor  pressure  (psio)  liquid  temperature  (°R) 

TBI  boiling  i urn j-erature  (°R)  pressure  (os ia) 

RH0  liquid  density  (lb/in3)  liquid  temperature  <°R) 

VIOCV  vapor  viscosity  (lb/ln-sec)  temperature  (°P.) 

CVAP  vapor  specific  heat  (Btu/ib-°P)  temperature  (°H) 

FKA  vapor  this  i  conductivity  temperature  (°R) 

(Stu/in-se*,-*  JF) 

CPL  liquid  specific  heat(Btu  r— VF)  liquid  temperature  (°R) 

liquid  temperature  (°R) 


ALM 


heat  or  vaporization (Btu/lb) 


Combustion  Specie*  Propertia# 

Subroutine  Property  (Unit*)  Independent  Variable  {Uhits} 


C0MP 

WTMINT 

TGINT 

VK1GAS 


DI7FU 


combustion  gas  composition 
(mole  fraction) 

combustion  pas  molecular 
weight 

flame  temperature  (°R) 


viscosity  (Ib/in-sec)  for  each 
species 

specific  heat  (Biu/moie~°F) 
for  each  spe  ries 

(thermal  conductivity  calculated 

frora  is  and  c  ) 

P 

molecular  weight  of  each  species 

calculates  diffusion  coefficient 
of  fuel  and  oxidizer  through 
combustion  mixture 


O/F  ratio  at  different  pressures 
Cpsla);  maximum  of  3  pressures 

O/F  rat?o  at  different  pressures 
(pais);  maxima  of  3  pressures 

O/F  ratio  r't  different  pressures 
(pels);  maximum  of  3  pressures 

temperature  (°R) 


temperature  (°R) 


Leonard  Jones  Parameters 

9-l  ,  atomic  radius  of  oxidizer  {hi 

« 

G„,  atomic  radius  of  fuel  (A) 

L 

Gj  ,  atomic  radius  for  each  species  (A) 

Cj/k,  interaction  energy  of  oxidizer  (°K) 

U/k,  interaction  energy  of  fuel  (^K) 

£ 

c,  '"k,  interaction  energy  for  each 
1  species-  (°K) 
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5.  INPUT  FORMAT 


PROGRAM 


PROPELLANT 

PROPERTY 

CARDS 


PROBLEM 

CARDS 


The  input  cards  may  be  divided  into  two  types:  (1)  propellant  property 
cards  and  (2)  problem  cards. 


The  propeU^nt  property  cards  contain  ail  the  data  shown  In  Section  3 
foi  the  particular  propellant  system  being  run.  The  properties  are  those  of 
tho  fuel  and  oxidizer  (both  liquid  and  vapor)  and  the  species  comprising  the 
combustion  gases. 

Input  Format  for  the  Propellant  Property  Caras; 


★Read  Propellant  Combination  Title  (A  Format)  in  Fields  1-10,  and  11-20. 

♦Read  NR,  NSP,  NTL(l),  NTL(2) ,  NTFLM(l),  NTFLM(2),  NT  use  110  format, 

♦Read  NSP  species  names,  (Fields  of  12),  6  fields  to  a  card,  in  same  order 
that  species  concentrations  will  be  read, 

♦Read  O/F  rati  '  (F  format),  seven  values  to  a  card,  until  NR  O/F  ratios  are 
listed . 

For  each  of  up  to  3  pressures; 


♦Read  pressure. 

♦lor  each  species,  read  compos  it  ‘on  as  function  of  O/F  (^even  values  per  card;. 
Start  new  card  for  each  species , 

♦Read  average  molecular  weight  versus  O/F. 

♦Read  gas  temperature  versus  O/F. 

GO  BACK  AND  START  NEXT  PRESSURE 


♦Read  in  gamma  versus  O/F  (at  any  selected  pressure). 

♦Read  NT  temperatures. 

♦Read  viscosity  (xlQ  }  versus  NT  for  each  species.  Start  ne 
species . 


card  for  each 


♦Read  specific  heat  versus  NT  for  each  species. 


Start  new  card  for  eaon  species. 


'i 


*mean 3  start  n«w  cord 


♦Read  in  molecular  waight  for  each  species.  Be  sure  to  maintain  correct 
species  order  in  all  input  data . 

*Read  species  interaction  energy  for  each  species. 

♦Read  propellant  molecule  Interaction  energy  for  ox,  then  fuel  (in  second  field) . 

♦Read  species  molecular  diameter  for  each  species. 

♦Read  propellant  molecular  diameter  for  ox,  then  fuel.  1 

♦Read  molecular  weight  of  ox,  than  fuel. 

♦Read  oxidizer  liquid  temperature,  vapor  pressure,  liquid  density,  liquid 
specific  heat,  and  heat  of  vaporization  (F  format).  Continue  on  new  card  up 
to  NTL(l)  valves  of  oxidizer  liquid  temperatures . 

♦Repeat  preceding  for  fuel  (NTL{2)  cards). 

♦Readgoxldizer  vapor  film  temperature,  vapor  viscosity  (xlO®)P  vapor  conductivity 
(xlCr),  and  vapor  specific  heat  for  NTFLM(l)  valves  of  oxidizer  film  temperature. 

♦Repeat  above  for  fuel  (NTFLM(2)  cards). 

♦Read  L00se:  110  format  f 

Fuel  decomposition  flame  exists,  L00SE  =  1  > 

No  fuel  decomposition  flame  exists,  L00SE  *  C 

Definitions 

NR  =  Number  of  O/F  ratios 

NSP  =  Numbar  of  species  >• 

NTL  =  Number  of  liquid  phase  temperature  | 

NTFLM  =  Number  of  gas  phase  (vapor  film)  temperatures 
NT  *  Number  of  species  gas  phase  temperatures 

(1)  ~  Oxidizer 

(2)  =  Fuel 

Follow  above  propellant  property  data  with  engine  data,  as  shown  on  coding 
sheet. 

Input  Format  for  Problem  Cards 
(See  chart  following). 


4  3 


Card  9  describes  the  geometry  of  a  conventionally  shaped  combustion  chamber 
according  to  the  following  definitions: 

Par  1  «  Z„,  inches 
c 

Par  2  -  er 8 ,  in8 

Par  3  «  (Par  4)  ZQ  +  rc,  Inches 

Par  4  «  (rc"rt^t’z<^ '  dimen8lonless 
Par  5  »  2w  (Par  4)3,  dimensionless 

par  6  *  2f(Par  3)(Par  4),  inches 


Input  for  -Sample  Problem 
See  following  pages. 
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IK PtJT  FORMAT  FOR  PROBLEM  CARDS 


■***  *J  t**+**3*f- 


INPUT  FOR  SAMPLE  PROBLEM 


^*0 

D2-14 

9  V 

10 

13 

7 

6 

4 

4* 

h20 

M2 

no 

439 

9  3 

0 

,5 

15. 

13. 

5, 

2, 

1 , 

.2 

.1 

'.00, 

,065 

,007 

,0? 

,coi 

,020 

,1M 

,408 

.560 

,61 2 
.113 

,a3l 

.161 

.27* 

.475 

,  JKR 

.219 

.090 

.046 

.351 

.033 

.360 

.47  3 

,i«l 

.  Jo? 

,36  9 

,348 

.340 

.  33) 

,005 

,0*5 

.013 

.533 

.479 

.33.J 

,0?9 

.002 

,011 

,020 

,090 

79,01 

28.37 

7ft,  5* 

22,7* 

10.  /5 

15.13 

i?.53 

U.62 

1903, 

11.3? 

7533  . 

3M66  . 

5290, 

527*. 

4018 , 

9693. 

U36, 

100. 

1  464. 

,004 

.015 

,01ft 

.179 

,410 

.560 

.612 

.113 

.  670 
.151 

.271 

.437 

,396 

,  <■  1  9 

.090 

,  0  *5 
.351 

.3^0 
.  363 

,  49J 

.  39ft 

.396 

.3*8 

,3*1 

.337 

.0*2 

,0*1 

,0i  1 

.  t>  33 

.*79 

,  3  3  J 

,  0  7ft 

.  302 

,cn 

.  0?t 

,007 

s 

>9,0? 

??,  i> 

?  6  .  >J 

22,90 

1  9  .  H* 

*5.11 

12,53 

’  1.63 
'  91  3, 

11.31 
?  i  3  3  . 

3  0  *>  5  . 

54VV  . 

6J9H  , 

*1  . 

?ft9*>, 

J 1 5  a . 

1  969. 

1  3  0  0, 

»**•«>•».-  , ..  ».i jtuwma Tm»*UKm&rtr.  ■•  ■■■♦-  .  \  *' 


.002 

,009 

,012 

,1*0 

.411 

.5^9 

,w,)9 

.526 

j  it  J 

.160 

.273  ' 

,447 

,402 

,??0 

.090 

.045 

.030 

•  391 

.360 

,3  73 

,394 

,308 

.370 

.347 

.348 

.337 

,037 

,007 

,516 

.477 

,333 

,076 

,003 

,005 

.010 

.023 

,004 

29,01 

28 . 3? 

26,6 

23,05 

18.96 

15.14 

12,54 

11,66 

11.38 

1993, 

2532, 

3672  , 

5558, 

5564, 

4014, 

2696. 

2145, 

1984  , 

1.310 

1.297 

1,253 

1,225 

1,23? 

1  .  748 

1.318 

1 ,  344 

1  ,3  53 

540, 

1080  , 

1V9J  , 

3060  , 

4860  . 

7020  , 

,419 

.659 

.976 

1,29? 

1.741 

2,20  7 

.50 

.79? 

1,179 

1,564 

2,108 

2.67? 

,614 

1.205 

?  ,  2  :)  2 

3,?5? 

4,693 

6.121 

,995 

l.:79 

2.367 

3,146 

4.24? 

5,37« 

i.to? 

1.731 

2,640 

3,511 

4,737 

6.00“ 

1,156 

1.91? 

2.85  3 

3,790 

5,1  30 

6,506 

.543 

5.173 

2.062 

2,932 

4,109 

5.283 

1,075 

1.795 

2,697 

3,571 

4,833 

6.13'I 

1,056 

1.747 

2,516 

3,466 

4,687 

5,944 

?,5 

? .  5 

2 , 5 

2,5 

?,K 

5.5 

7,0 

7,0  5 

7.35 

7,9-9 

8,62 

9.11 

5,08 

4,73 

10,22 

11,75 

13.? 

1  4. Cl 

7,0 

7  .  ?  4 

9,0 

9.51 

9  ,  *5 

«  ,  90 

7.0 

7.04 

7.45 

8 .09 

8,73 

9.1 

7,07 

7.7 

9,47 

8,89 

9,44 

0 . 9? 

5,55 

10.73 

14,3 

10,53 

18.43 

'9.36 

7,17 

7.5 

9,33 

8 ,59 

8.97 

9.12 

S>6 

5.09 

5,3  2 

6.01 

6,0? 

8.12 

1  , 

3. 

16. 

28, 

1  7. 

<2  . 

1  7  , 

30, 

16. 

37, 

59.  / 

909. 

71,4 

7  9,8 

10  6,7 

668, 

116,7 

106.7 

331,9 

502,9 

2,71 

3.83 

2,64 

3,8  0 

'  ,  1  5 

7,4’ 

?.9 

3,49 

3.05 

3,  969 

4  .  C2? 

*6, 

7  ?  , 

40  0  . 

.27 

i .  5  j 

,35 

4  35  , 

440  . 

,  9  5 

.  C  5  5  n 

.  <5 

A  t 

• 

4  40  . 

'  .  94 

.  :  5  4 1 

.35  3 

4  2?  . 

c-  2  *•  . 

‘1.5 

. r  9  <  y 

,372 

*12  . 

47 


969. 

SO. 

.09385 

,  J6fl 

403  . 

609. 

60. 

.  0497 

.414 

39-9, 

*60 , 

100. 

.0463 

,4  M 

370, 

689. 

330  . 

,04  Vi 

.543 

365. 

729. 

660. 

,0  \H'J 

,70 

34  8  , 

760, 

1  200. 

.  0324 

,87 

316, 

490. 

.091 

,035 

,  6*6 

706, 

509- 

.GO 

.  0  366 

,70 

680, 

*59, 

.46 

.0137 

.72 

650  , 

*00. 

1.9 

*  0  ’  48 

.74 

62 n  , 

690, 

6.3 

,0338 

.75 

59n  , 

766, 

19. > 

.0328 

,77 

6>n , 

790, 

37. 

.C3i7 

.76 

510, 

*09. 

75. 

.0316 

,80 

470  , 

"40, 

160. 

,0792 

.82 

430  . 

900. 

250  • 

.0341 

,83 

390  , 

990. 

430. 

03  70 

.85 

37n , 

1090. 

660. 

,0360 

,«7 

340, 

1090, 

960. 

.0??0 

,89 

ill, 

400  , 

,6?4 

,0  94 

,231 

600, 

.92 

.234 

,238 

1000. 

1  .57 

,44j 

,253 

1600. 

2.34 

.1,76 

.269 

7400, 

3.21 

1  .49 

.286 

5200, 

3.91 

.91 

,297 

4000, 

4.57 

2,4/, 

.  300 

400. 

.4? 

,14 

,-<2 

*00. 

.80 

,50 

1200. 

1.19 

.49 

,61 

?ooo. 

1 .95 

.  9-> 

74 

2899. 

<r.71 

1.3 

.02 

4000  . 

2.83 

2,4 

.90 

N244» 

i 

8**h  st*w;uty 

/  JM69  HuN 

183 

f  DOU*U  T*N?m4  )  I  NjCJ^- dO 

"EUV« 

vi\j  93  r: 

OK^h'T 10MS 

IC4P6 

5TM 

0,9 

.05 

.0331 

.01 

1.6  "  ,  0 

0 

«  1 

1  1 

Hi. 

.2 

4  . 

12,46 

128.6  1.4 

12.99 

646, 

5j4  , 

.9H565 

7,3  5 

11,24 

1  60  9. 

5J3  . 

,0Oi»RO 

3.7  5 

128.4 

48 


PROGRAM  USTING 


6 . 


PROGRAM  ST=At)Y<  J  PUT, OUTPUT,  TaPE5«!NPUT,TAP(EA«0UTPUT> 

SJRFTC  TRSl  li$t,bef 

000003  COMMON  /S/rC»NM/R/8D/C/TA/FP/AP, AU.BCD*#RCDY.lOft 

000003  COMMON  /LlYKA/  LOOSE 

000003  CO^hJN  /PR3S/  PRPSYSM)  /SPEC/  NSP,  SPNA*E(2.  15  > 

003003  O!«£YS!0N  rC(80 ) , N*(20 ) • BD( 1 0 , 20 . 1? > . T * ( 1000 ) . ABI ISO > . 

1  AJ(iOO,90>.8CD*(7),BCOYC3*,LOP(9).T!TLE<24» 

?,£H(10,20*?)»  OH<10,20#2), YLA(10.20,2),  DTdn»?0«2) 

3,Vd0.2D,2>.  OV  C 10 « 20  <  2 ) ,  *T*0L(2> 

000003  FOU! VALENCE  <FC(14),  TSTSP),  (N4(5),  RPTN ) •  (BRd,  l,  1>,  c*d,  1, 

11)),  (R3U.  1.  3 ) ,  D*(l,  i,  i>>»  (ROd.  1#  5>,TLAM#  1,  in#  tRDll 
2,  1,  7),  OT(i,  1,  1)),  (80(1,  1*  9),  V(l,  1,  ID,  (80(1,  1,  11>,  0 

3y(i*  l.  in,  <rc<i),  n,  (nn(15),  mplot),  <t*{4 tu ,  utmod 

000003  COIHON  /TRFGAH/  Ff»A*4T  d3 )  /OF/  R(l5>.  N»  /CMPTAP/ 

1YS°( 15,  15,  3)  /PRSSR/  PRSS  <  3 )  /PVY/  PV(i5,  ?)  /VI Xf  Tt(l5,  ?)  /DT 
2L/  NTL(2>  /TFL4X/  TFlLUlS,  2)  /YJSCVY/  VVIS(l5,  2)  /DTFL*/  *TFL*M 
3?)  /RMOr/  9HOLd5,  2)  /  K  V  A  P  Y  /  KYAP(l5,  2)  /CPVPY/  tPVAP(l5,  2>  /CP 
AY/  C » ( 15.  2)  \kny/  LAH(i5,  2)  /WTH/  w*T*Bd5.  3>  /TGI/  TC1I15,  3) 
5/TA8'/  f  T  A  3 f  ”  5 )  /V*<SGS/  VS»(15,  15),  CPT*3d5,  15)  /0T/  NT  /MFFUS 
6/  EP3I(15,2),  S 1  SM«  ( 15 ,  2>,  FPS0F(2),  SIGOH?)  /■**/  WH(i5) 


000003 

9£AL  K»UP.  LAH 

000003 

B6*0  (5,  601)  PPPSYS 

000011 

601 

r0BM*T  ( 4 A  6  ) 

0  0  0  311 

»E*0  (5,  60)  NR.  NSP,  VTL,  NTFl*),  VT 

000027 

60 

FORMAT  ( 7 !  10  ) 

000C2? 

P£»D  (5,  603)  (SPNAHEd,  I),  SP*ANE(2,  !),  I  *  1,  \5P> 

00004* 

603 

FORMAT  (12*6) 

000046 

ofe*0(5.JQl><  P'D  ,  1*1, NR  ) 

00 3 061 

no  15  I  NO* 1 , 3 

00  03  "' 

b£  *  U ( 5 . JO  1 )  PR5S(^D) 

00  00  7 

"0  1?  1*1, YSP 

00007? 

P£  *  D ( 6  « J01)  (  Y?P<  J, I , J  NO ) ,  J*1 , NR  ) 

G  0  0 1 0  v 

12 

C0nT!nJ= 

000  1  1'’ 

t-E*0(5.J015  (  wmT*m(  j  ,  i no  ,  !*!,*;<  ) 

0  3  0  1  25 

9  E  *  o  ( 5 ,  ?  o  i  >  (  rru  * ,  i  nd  >  ,i*  i,NB) 

300141 

15 

CO'iTiNJr  , 

000143 

BE  4  L ( 5 , JO  1 )  (  FG*^Y(J), J*1,N=) 

030156 

BE  13(5,  JO  1M  .  T  A  R  <  !  )  ,  !*l,vM 

0001 71 

PO  17  1*1, NS® 

3  C  0  i  7  3 

fit  *  D  (5,901)  (  VSP( j, I ) ,L«1 ,NT! 

C  0  0  2  0  ( 

17 

CC'dlN.'E 

000211 

no  i?  :*i,ys» 

03021" 

BE  » 0  (5,901)  (  CsTtfi( J,  -  i,  J*1,ST) 

C  ?  3  2  2 

19 

c  0  ^  T  i  N  J  = 

030233 

bE*C(5.joi)  ;  -»<?),  i*i, sop) 

050242 

B;  *D  !  5 , JC1 )  (  EPS] (1 , U . !  -1 .NS°5 

000255 

BE  *  D ( 5 • JO  1 >  EPSCF 

03 0263 

s  t  4  U ( 5  » J01)  <  S:Gm*(:,15,  fxl.sS3) 

0  3  0:?'- 

B».  *D<5,JC1)  5  1  G  0  r 

3  0  0  3  3  « 

sF.*0(5,J0i)  »tnqu 

•  '**<*M*W'*w*r't 


06131* 

•••314 

mm 

imo 

01333ft 

606330 

300334 

000401 

030403 

000411 

Q004i< 

060411 

000411 

000411 

03Q41? 

000413 

030«a* 

000416 

000420 

30042? 

090426 

000431 

0*04J? 

£0043' 


000434 

00043? 

003440 

000441 

06044? 

000443 

00044* 

000446 


no  20  1*1*2 

saCT«3Tl<  J  ) 

*6*0(5, »02>  <  TLIJ.P  *  PV< J. P *»H0L( J. I ) *CP( J. 1 > ,  LAN(J.H« 

1  ) 

20  CO'iTiNJ; 

00  23  1-1.2 

uTfCT-  *TF.MI) 

0140(5.003)  (TMIKIJ.I  ) .  W 1 S  <  J.  I  i,*VAP(  J,I>,!*PVaP(  j,  !  ) , J-l , NtFCT > 
23  C0*T!*j= 

*6*0(5,904)  tOOSt 

901  ro»*4T(  7110*4) 

902  rQ**lU  5610,4) 

903  F0**4T<  4610,41 

904  r09H4T(|lO) 

5  CALL  MiiO 

CALL  39G 
CALL  *$ET 

to  CALL  *«Oa(NSP> 

CALL  OifUVfNS*) 

call  v*D*<D*»§M,n 

call  *A04(f)T,Ti a, 2) 

CALL  4*04 ( 0V , V  ,  3) 
call  sum 

IF(4>tN5  40.36.4S 
36  C*LL  MINT 

C  CENTS'). S  STEP  PPJNT  4*0  S?0«ES  PLOT  V*luE> 

r  m6*T  T=ST  “00  T I -6  STOP 

I F  (  T  -  TS?0*»)  40,  39,  3ft 
3ft  t f ( NaL3T  ,4E,  0)  GO  TO  5 
CALL  P.3TT 
C.o  TO  5 

40  c  ALo  MITT 
50  To  10 

12346  CALL  fc<|T 

r  NO 


SO 


function  Aon 

®88®®S  COMMON  /  sDUPHM  /  PA«1#  PAR2,  P*P3»  PAR4*  Pfc»5a  PA*« 

000003  JF  (X  ,3T,  PARI)  so  TO  I 

00000*  a  ■  par 2 

00000?  RETURN 

000007  1  A  8  3,141*  •  (PAR3  -  PAR4  *  X)  *•  2 

090013  RETURN 

000014  end 


FUNCTION  AP  (  X ) 

GOU0O3  COMMON  /  ADO*«M  /  PAf<l#  PAP2*  PAP3,  PAR4,  PARS,  PAP6 

O0GCO3  IF  it  ,  3? ,  PA&l)  QO  TO  1 

000006  AP®0 « 

000806  RETURN 

00000*  1  AP  *  P4R5  •  X  *  Pa«6 

000912  RETURN 

090013  END 


090005 


00000*5 


090005 

00PCC5 

090011 

000021 

000021 

000022 

000043 

000050 

00006? 

090063 

00010H 

000110 

00011? 

000115 

000120 

000123 

000137 

000143 

000151 

000164 

000170 

000179 


SUHROJT1N?  BUlLtJiI) 

D? PENSION  Fg(80)jNN<29)»9n(10«20#12)»TA(1000 >.BFTa<?0,2)«B5<29,?) 

lfcFNUM(20.2).UW(20#2),DM(18.28,?>,DTH0.20,?>.V(10,?0.2). 

2  TLUO.20.21  ,EMUu. 20*21.  FNUMI20.2) 

EQUIVALENCE  CT4(J0),P), <TA(39?) ♦RV),(TA(474),FKAV). 

1  ( TA ( 26 ) ,  TWOP I  i »  <Ta(473!.CPA), <T4<482),QJR), (TA<2A4).TG), 

2  (TA(64?),5>E7A),  (TA(  96), RS>,  C T A C4J3 ) , FNUN ) ,  CBDI401  >,DM>, 

3  iCD(l?Jt),D?>.  (BBC801),  TU),  <8DU),EM>,  C«*DU601),V) 

4,  ( T  A ( 4  52 ) •  G JR )  , (TA(393),FNUH)  ,  t T A < 216 ) , Ww > 

COMMON/E/pc,NM/B/BO/C/TA 
150  TR  «  73UM) 

wRITg( 6,2000)  J,  I 

2009  roRMAT  (IK,  SHOORQP,  12,  IS.  10«IS  BOILING) 

<1  ■  3 

wl  ■  A4!N1(ABS<-DM<2,  J,  I)  •  VCl,  J,  m,  l.E-6) 

RRS  ■  RV*RS( J, I )**2 

COE  *  e‘NUM<J,I)*FIKAV*RS(J,  D*TrOPI/CPA 
COE  «  ABS(COE) 

152  R  *  CO?  *  ALOGU,  ♦  CPA  •  (TG  -  Tfl)  /  (AL*CTB.  I)  *  GJR  •  Iwl  /  RR 
IS)  •«  2)) 

WR  a  AR3<(W  -  wl)  /  w) 

Hi  «  <1  *  1 

IF  (41  ,GT,  15)  GO  TO  154 
133  wl  «  (4vwi>/2,0 

|F  U'R  , GT ,  ,02)  GO  TO  152 
154  DM  1 1 , J, I )  *  -w/VU.J.I) 

OTci.j.n  *  0,0 

TL( 1, J*  t )  *  YR 

mETA(j*  I  )-  *  w/«S(j,  I  )/FNUM(  j,  I  j 
w«( j, n  *  w 
RETURN 
EM) 


51 


COCCI* 

0  0  0  6  i  4 

70 

tflRi-l  .5-4)30,30, 40 

00002* 

3" 

CD  *0, ?71* ^fc** , 21 7 

0  30025 

9ETy?M 

000029 

40 

CD*  2 . 

000027 

RfcTUPN 

000027 

ENH 

0  0  0  0  0  5 
0  0030? 

ooooo? 


00000? 


00000? 
000005 
0006*7 
000011 
000013 
000014 
000020 
000025 
0  000^9 
00004? 
000057 
000061 
00006? 
000066 
000076 
000107 
000120 
0001 J4 
000147 
00015? 
000161 
000167 
000173 
52 


SJ3HUJT1NE  nfcHlV(MSP) 

COH^OM  /6/f C, VM/B/B0/C/T4 
COMMON  /L 1  \IK A/  LOOSE 

nJMtSSION  rc(B0).MH(20).  PDd".  2'1»  1?>.  TAd000  ).NSET<?) 

1,  6NMP(2>»  KniSC(2o#2)i  Vl)(20#?)»  E  M  ( 10  ,  20  *  2  > .  T|.  ( 10  ,  20  .  2 ) 

2, V<10.?3.2>.  M  S { 2  0  «  2  )  i  WTMPL«2»»  PR<?0  .2).  SC(20.?).  RE(20,?) 

3 ,  REP(  20,2),  FNJH(2U,2),  FNHH  ( 20 , 2 ) ,  !)*<<  1  0 . 20 . 2 ) .  nT(10,?0.2) 

4 , DV ( 10  *  20 ♦ 2 )  #EM!NT(20.2) 

PQU1VA  =NC=  (TA(461.NSET>.  ( T A ( 4 79 ) , E*MP > ,  (  T  A  (  27  fl )  ,  K  P  !*C  ! 
r,(TA(257>.4AT),CTA(30),P).<TA(476);  PA).  (TA(475).  CPA) 

2,  (  TA ( 392 ) i  «V)  .  (TA{?11>.  VD>.  ( T A ( 2B) . F TP  1) . ( T A ( 29 ) .  TH90> 

3,  (90(1),  E4)  »  (BP( ► 0 1 ) .  TL).  ( HQ ( 1601 ) .  V)  ,(TA{96)»  RS) 

4 ,  (  T  A  { 2  5 )  *  3),  (TA{48),  WT  Mf>|. ) ,  (  T  A  (  47*  ) .  U )  .  (TA(391>,  9H0VS) 

5,  ( TA( 443) . 4TML ) •  (Ta(474),  FKAV),  (TA(473)»  nPOT > , ( Y  A ( 48  7  ) ,  DR» 

6 »  <  T A ( 5 ?7 ) #  SC  )  #  ( T  A ( 567 ) i  96).  ( T A < 407 ) . REP > .  ( Y A ( 393 > . TM"H ) 

7. <TA<433).-NUM).<BP<*Ol>.OM),  1BO(120D.nT)  .(911(  2001),  OV) 

8 ,  ( T  A  <  4  9  4  ) ,  RH0G),(7A(264),TG>  . ( T A ( 56 ) . Em I  NT ) 


P2  *  P* 2 . 

120  no  200  1*1.2 

N»NSFT(  !  ) 

ENMP(  1  )  *0  . 
no  200  J*1,'N 
Knisc<J.n  *  o 

ir(SH(l»J*  1))  200,  200 i  130 
130  T9R  *  (  TG  *  TLd,  J»  l  )  )/2. 

CALu  VKSGAS  (RAT.P.TBW.VjSB.rKB.CPB.NSP) 

PA«A?1V1( ,99*P,PVAP(TL(1, J. !>»?>) 

PA2P*P A/P2 
DP *2°*1  ,-P*2P 
CPA«CV**( T3H,  I ) 

9H  *  R40<TL(1,J,I)  »!) 

9V  «  WTM0L( 11*P/R/TL* i. J*  l> 

VU( J, I )  *  U  •  V(1,J,1 ) 

PS ( J. I )  *  (EM(l.J.I) /FTP !/PH)**THRP 
RrCJ.  1)  •  9S(J.  I)  *  PHOVS  *  A9S(VH{J,  !)) 

RTH4V  *  PA2P*WTM0l  (  I  >  ♦  QP(  p-kITML 
VISA!/  *  V  JSCV  (  TRR ,  I )  *PA?P  ♦  ->A?P*V!SB 
F<AO  *  3A?P*FKA(TflR. !  )  ♦  |)PA2P*FK9 

mo  a  v  *  p*wtmav/r/tbr 

C3AV  ■  (PA2P*«TM0L( I )*CPA  ♦  DPA?P*WTML *CP9 ) /WTM A V 


■> :  k  ? 
:  o  o  ?  i  * 


C  •0  0226 
000231 

0  0  0  ?  3  7 
000244 
000261 

973 

000277 

797 

000277 

079 

000305 

000321 

000333 

000335 

000343 

r 

160 

000343 

000354 

000356 

000365 

161 

Q00367 

162 

000375 

000400 

163 

000402 

000404 

149 

000406 

000411 

150 

000413 

000415 

155 

000417 

156 

000451 

200 

000457 

000457 

OwGAV  •  *951  VOt  J,  I)  ) 


:  ;  *  c  ,  *  *  S  i  j  »  n 

i*‘p<  X  **  A  /  2  , 

tius.  yi-ru  craTiP.tpr.paip, i .nru.NS**) 
D»OT  «  3Fu#P/TBP 

®H ( J  » 1  >  «  CP*V*V1SAV/FKAV 

SC(J.n  •V!S*V/RMOAV/OFU 
JF  <P9(J.  I)  ,  GT ,  Oi  ,  AND .  SCCJ,  I)  ir*T, 
WRITE  (6,797)  PRC J, !>,SC< J» D 


0.  ) 


FO^MATCHO.ISHSTATEMFNT  97B  PR  *,  F10.3.6H 

RTRE  «  SORT(REP(J,  !>> 

rNjM<J.l)  «  2.0  ♦  ,**SCC  Ji  I )**?hrr*rtre 

FnjMC  J,  I  )  «  2,0  *  ,6*PRC  J»  !)**TMRn*RT«E 

CmEC<  FOR  ENTRANCE  TO  DISSOCIATION  LOOP 
!F  (LOOSE  ,EO.  0)  GO  TO  149 
GO  TO  (149,  160),  ! 

CONTINJE 


GO 

SC 


9EE  *  2,*  RS(J«l)/(FNUH(J«n  *  ?.> 

CALL  U/CHKCKNOW) 

00  TO  (161,  162),  KNOW 
9EE  *  1.EJ0 

RL  «  RTC  RS< J,  !  ) ) 

[ F ( B  =  =-  RL)  149,163,  163 
CAL.  T4FL  <J,P 
30  TO  1*6 
PMOP«0 • ?9*3 

IF  ( B  A  ,LT,  PHOP)  GO  TO  155 

cal.  boil  fj»!> 

GO  TO  156 
CALL  (CAPOR<J,I) 

nv(l.  J,  I)  «  RHOG  *  ,375  /  RSU*  I)  *  ARS  ( VO  ( J, 
1*  VD(J,  I )  /  RH  *  CDCREfJ,  I)) 

CONTINUE 

Rfc  T  JRN 

END  _ 


000003 

000003 

000007 

000011 


FUNCTION  FGAM(OFR) 

COMMON  /OF/  R ( 15 ) i  NR  /TRFOAM/  FGAMT ( 15 ) 
CALL  LN3RIN  (OFR,  FGA  W ,  R ,  FGAMT,  DFDT,  NR) 
RETURN 
END 


/  V  5  5  A  V 

TO  979 
* , FI C , 3  > 


I  ) )  /  V(l.  J,  I 
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0  5  C  5  1 1 

1  sifc  ung*?n  u#n,»ft*G,x  .y,oyo*iN) 
ni -6vSiON  * ( i ) , r ( ii 

o  e  o  o  1 1 

r 

MM  *  M 

090011 

r 

jf(mm,m  =  ,i>  r, o  to  l 

O00014 

YAH3*Y(1) 

000014 

D YQX ■ 0  « 

0  9  0  0  IS 

ft* 

HfcT'JHM 

0  0  0  0 1  ft 

i. 

i 

XVSA'JG 

000017 

1  r  (  x  <  M  M  »  ,  L  T  ,  X(D)  GO  TO  3 

0  9002.1 

no  2  j*2*nn 

"30024 

1  e-  <  x  <  j),G=.xv>  GO  TO  ft 

000027 

2 

rnviT  i  n  j  = 

0000 32 

GO  TO  5 

00003? 

3 

no  4  J«2#NM 

000034 

ircx(J).Ls.XV)  GO  TO  6 

000037 

4 

COMTINJ5 

00004? 

5 

J«MN 

000044 

6 

I  *  J-l 

00004ft 

nvox«< yc J)-y< i ) )/( * ( ji-x(  i  n 

000054 

ya9g»y( 1 5*DYnx*<xv*x<  1 1  > 

000061 

c 

e?£TU4N 

000061 

FNH 

54 


r 

c 

000007 

000007 

mono 

00001*5 
000017 
000021 
000025 
00002*5 
00C327 
000021 
000034 
000035 
000040 
000043 
000044 
000046 
00005? 
000060 
00007? 
000076 
000103 
0001  06 
000110 
000111 
000114 
000117 
000124 
000126 


j-.t 

fj-PjTts  R*bli  Of  N  DROPS  *B«U7  which  l/u  Th  OP  « «F  h*5S  ;  5 

I  v  ’*?  LOG  AO  \  THm  J  CO-  NORMAL  niSTRIBUTJON  wtTw  *US*-W«Um 

ra-ijs  r*  tur  standard  deviation  in  sig,  *p.  the  inverse 

PL'lb*rj!.lTy  rjvCTJON  IS  APPROXIMATED  BY  A  RATIONAL  f  R  A  C  7  J  i  •  * '  • 
HIMfeNSlDN  RS<20) 

OfcL  st,/  PLOAHN) 

S!3lV  *  AlDbfSIG) 

3*-D?L/2, 

P  «  s  v  ♦  4 

ME-P  *  <  N*1 »/?  ♦  4 
J  *  MP4 

no  4  I  -  i , NP4 
iru-2)l0»15,ll 

10  PtO*  0 i 5*D:L 
GO  TO  2 

11  ird-5)  15,15*20 
is  Ptp*o , ?*rifc^ 

GO  TO  ? 

20  psP^DEi. 

30  IF  (1  , jT ,  MEO)  GO  TO  5 
?  t  s  sort ( - *log f p  *  p>> 

Xp»t  *12. 30753*.  27061*T1/U.  *T*<  , 99229* , C 44Bl*T >  > 
esiGX  «  E  X  3 ( S I GLN  *  XP) 
f  F  < 1-314,3,40 
40  IF< 1-534,49.50 
49  p  «  06./2.0 
GO  TO  4 

5C  J»N  -  T  »  9 

3  RS< J)*R4*ESI5X 

4  RS( I )»94/feSlGX 

5  »ETU9N 

END 


000011 

000011 

00001? 

000014 

000016 

c:::2o 

000021 

000024 

000031 

000033 

000036 

000043 

000044 


Subroutine  m asset <e mo, rho.rs.emdot^.opopn) 
n  I  yE  MS  I  3N  =liO(20)|RS(20),DROPN(20> 

GRwO=4, l«B790204*RHO 
FR*C  *  =MD3T  /  FLO»T«Nl 

TRACI  *  fmac/6,0 

N)P<  s  \i*4 

r.u  2  \  *1,1 

FHD( I )  *  CRHO*RS( I  1 • * 3 
5  DROPM( I )  *  FRACl/gMOl  I  > 
rn  5  i*5, N34 
CMD< i >  *;RHD*RS(  I  j**3 
5  nw'jPMt  I  )*f  RAC/EMOC  I  ) 

return 

END 
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0  0  0  3  0  6 

0  0000  6 
o  e  i)  o  o  a 

SUBROUTINE  MUM  YP,  V  ,K«  ) 

DIMENSION  >  K  210  ) ,  n’lOl,  FC<80),  Nm<20),  u  T  (  3  ) ,  NST  1  2  ) ,  NE  D  (  ?  ) 
COMMON  IE/  FC.Nh 

EQUIVALENCE  <Fr<3),H?;  (  F  C  <1 1 1  ,  W  T  ) ,  f  r  r  <  1  C  >  ,  El),  <FC(9)t  FmAjO, 

1  (FC(8),  feXIN),  (  NM  (  2  )  <  t  NPR  )  ,  ( NM ( 7  )  »  NSM),  <NM<4).  M  ^  A “  )  # 

2  { NM ( 3 ) <  malR)  , <NM< 17)  ,NSt ) , (NM( 19)  ,NEn> 

r 

r 

MtiAMt-1,0,1  INDICATES  PREDICTOP-PUNGE-^UTTA  OR  cnRRfcCTORDHASF 

c 

HALP'INDICATES  PERT  OF  RUNfifc  KUTTt  PHASE 

r 

HSP-  PRINT  SWITCH  FQR  ERROR  INDICATION 

c 

FI  -  C I  STAINS  MAXIMUM  ERROR  FOR  EACH  CVClE 

c 

VP  *  address  of  derivative  array 

c 

Y  *  *DDRESS  OF  The  ORDINATE  array 

0  9  0  0  []  (■ 

c 

IF  ( mg  A 1 )  40,  10,  60 

r  this  ro jt i ne  emplo' s  a  fourth  ordeo  pungf-kutta  st*rt?& 


r 

000010 

10 

IF  (mai.3  -  2)  20,  15  >  15 

0  0  0  0  13 

(X 

K  *  4  •  M  A  i.  P 

080013 

DD  15  K.  *  1,2 

nnnoiT 

LS  ■  NST(KL) 

000021 

LE  ■  f,6D  <  XL  ) 

000023 

00  13  l  •  LS , LE  » IT 

u00025 

j  *  1  ♦  K 

000027 

Vi  1*2)  «  Y( J)  *  H*  y P i  I  ) 

000034 

IS 

CONTINUE 

000040 

CO  TO  99 

0C014S 

20 

no  22  «L  =  1,2 

000043 

L S  *  NST(KL) 

000043 

LE  *  NED(kl) 

000047 

00  2?  i  S  LS.LEilO 

000051 

Y ( ) ♦ 9  )  *  v  <  1*3)  ♦  H/6 , * ( YP ( I) *2 , * ' VP l  T *1  )  *YP ( I *2 )  »  ♦  YP(I®J!) 

000064 

22 

continje 

J0007P 

r 

GO  TO  99 

c  4DAHS*-<tSHro»TH  fjfTM  ORnEP  PPFOiC^OP  FORMULA 


i-  *«»***»  t  »*»  *  »* 


''10070 

40  on  42  <L  *  1 » 2 

U  0  C  0 7? 

L?  *  NST(XL) 

0  0  0  0  7  » 

LE  *  NED (XL) 

003376 

DO  42  I  «  LS , LE , 1 0 

UDUUIO 

T  c  I  ♦  9  )  *  TU)  ♦  H*(?,64027777777*YP!  I  )  *3 , 85277  7  7  7  7  7  7  *  Y  =  I  1  ♦  1  ) 

1  ♦  3,53333333333*YPn*?>  »  l , 76944444 444. Y6 < t *3 ) 

2  *  ,34<illimil*YF(  1*4)} 

000120 

42  CCNTINJ: 

000124 

50  TD  99 

f  *UAHS  HTUlfCH  flFTM  ORDER  CORRECTOR  FORMULA 
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O  0  Q 1 2? 
0  00127 
0  5  0  V  3  i 
0  uti  1 33 
00013? 


000154 
C  0  0 162 
000164 
000170 
000173 


Q0017* 
090200 
000201 
0  0  C  ?0  7 
300207 


000212 
000214 
09021? 
0  0  0  2  2  3 


000223 

00023? 

000236 


r 

c 

c 

m 

c 


60  M  n.  *  1*2 
1 5  *  NST<«l> 

Lc  •  "te'll  *L  J 

uJ  93  1  *  lSiLEiIO  „  „  .  , . 

r<n  «  7(1*1  >*HM  ,348611 1!  1111*VP(  J)*  .8972?222?222*VR( I *1 > * 

4  n  t*.6?666?6?*VP(  !  *2  >  *  ,  14722222272*YP(  1*3) 

2  -  '  ,  0263BS98B»888*VrM*4n 

6?  ASS'  U  D  -  Y(  1*9)  )*HT(KK) 

I r ( v «  I  )  )  70  ,  80 i  70 

79  fc9  «  =9/  ABS(  YU)  > 

80  1 F { ER  ♦  CHAX)  65*  95,  9? 

9?  If<ER-  EMIN)  98,  87,  87  ^  ,.,*«•*#*** 

RfeL A  T  I  V  ?  ERROR  CW£CK*fe)RANCW  TO  99  INDICATES  FRRQR  SMALLt^  TMAN 
ALLOXAS~F  5R90R»ADniMQ  ONE  TO  INDR  INDICATES  VARIABLE  JlThilN 
ERROR  A. LOWED  ^  **#•••###** 

f.  f  #  £  i  *£&**©#♦♦#•*****•****•*•***  ****** 

9  7  I  NOR  «  !NuR  ♦  1 

ir (NSW)  sa.  *a,  «8 

5 ft  wRITE  <  6.2',  XX 

2  rORMiT  il‘!3.  8MVARURLE,  !3,  13MWITHIN  LIMITS) 

GO  TO  98  # 

ONE  HUNDRED  IS  SUBTRACTED  7 C*R  EACH  VARIABLE  LARGER  THAN  The 
ERROR  >_ !  m  I  T S  ^  ***«•*#*••* 

55  INIJR  «  I NDR  *  100 
1  F ( NSW  )  97,  98,  97 
97  wR  I  TE (  6,3)  XX 

j  FORMAT ( 1H0 ,  IHVARJABLE.  ! 4 ,  20MOUTSJDE  ERROR  L ’ M ! T5  ) 

Pi  CONTAINS  MAXIMUM  ERROR  OCCUR  I NG  DDRtNG  Tw  =  CYCLc 

.*«* i*  ****.*♦***»«*«*****  ******* »•♦**••****••*****•* *»•**•*****« 

9*  El  s  a  wax l { ER , El ) 

99  RET  JR  N 

END  _ 
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R  „  a  R  "5  U  t  I  \  e  '  9  J 

0  0  0002  D1*ENSI  3*  rC'#Ci,  1  “  ■  2 f  )  *  !  J  ’  ,  *  S 1  T '  2  ■'  i  *. S  ’  '  .5  :  ,  s 1 "  (  ?  !  i  ? *  •  i  “  * "  ; 

n n o o o ?  coh-om  /]t>  fc,  /c/t* 

38000?  £  C'J  I  V _E  N  „  £  (  F  C  ( 1  >  *  T  >  4  !  r  r  (  ?  )  ,  » “  ’  i  ,  t r  C  <  J  •>  .  -  )  ,  <Fm<4),hC), 

i  (fC<5),  hhin),  t  ret  a ! ,  mm*w  > ,  ( rr  <  ?  ; .  *ZD' ) .  t r  m  is ) ,  £  * i  v ) , 

2<ec<9>,£max),  irciiii.m,  tret  ie )  ,<-i  > , 

3  (MH(1»,!K).  :mm(2).!n0»>.  (  WM(  3 )  ,  M*L e )  ,  (  nm ( 4 ) , “GAM ) . 

4  <MH(b)|HPTM,  <  MM  (  6  1  .“PTS)  ,  (  NM  <  7  )  ,  A'S*  >  .  <  Nm  (  9  )  ,  NCOU  )  , 

5  (MM(9S,MP),  (  NM  (  i  0  )  i  N  V  )  ,  (NM(ll),  MO) 

6,  ( TA t 46) , NSET : ,  ( N* ( 1 7 ) , MS T  ) ,  (NH(IQ),  NCR) 

e  »****«.****t»*.»MM«M*«.#...».»*.*.*.*.»**.**#t.*tf.i..«.»*. 

r  DESCRIPTION  of  The  listed  v  aw  t  »rle$ 

r  T  •  THIS  CELL  CONTAINS  CURRENT  INTEGRATION  TIME 

r  RKT  -  START  TIME  OR  PRFVIUUS  BEGtNING  OF  RK  TJHJ 

r  H  *  CURRENTLY  used  ST£S>  site  IN  COMPUTING 

c  HO  *  STORED  STEP  S I  ?E 

0  HZ02  *  HAL r  Of  STORED  STEP  SIZF 

C  HHIN  -  41NIMUM  STEP  SIZE 

r  hmax  maximum  allowable  step 

C  FMIN  -  EM A  X  M 1 N  AND  MAX  AOlOWARLE  FRRQR 

r  w  -  array  of  WE  I  wn'i  -»  .0  WEIGHT  ERROR  CINSIDERATION 

c  I M  -  NO  OF  GOOD  POINTS  FROM  R,K  START 

C  INOR  »  INDICATOR  FOR  ERROR  OUTSIDE  OR  wJTmIn  M t N  max  TOLERANCE 

r  MALP  -  COUNTER  for  R,K  INTPRMEDIATF  POINTS 

C  hG A M  -  PHASE  INDICATOR  *t»PRFU!CTOB0»R}K  1,  jORRECTOR 

r  MPTN  -  PRINT  COUNTER,  CURRENT 

r  hp T S  *>  TOTAL  no  of  POINTS  in  PRINT  INTERVAL 

C  MS-  -  ’HINT  INDICATOR  IN  NaDM  ROUTINE 

r  NCOU  *  TOTAL  NO  OF  COMPUTED  POINTS  DURING  INTEGRATION  CYCLF 

C  hp  -  POWER  OF  2  variation  FROM  MMJN  TO  hma* 

C  HMAX,  W»,  NO,  NSW,  NV  AND  W (  T  )  MUST  EITHER  BE  READ  INTO  CORE  0° 

C  I N I T I  A | Z  E  0  BT  AN  ADDITIONAL  ROUTINE 

r 

00000?  HMIN  «  HMAr/?,**Hp 

0  3  C  3  0  7  HO  »  H h  i  N 

080010  hZ02  *  W0/2.C 

00001?  h  «  HZDZ 

OOCOIS  PwT  «  T 

0  0  0  0  1  -•  £1*0.0 

r  •ft#**ttft«*«*»»*«t«**»t«t«#«**i«**»*****»*6***«»«*******t#***t«« 

r  FIXED  POINT  INITIALIZATIONS 


r 


ocoois 

1-  *  0 

0  3  0  01* 

HALP  *  4 

0  u  *  *  * 

'  ’AH  t  0 

OCwO 20 

m=»Tn  * 

0 

000021 

HP  T  S  « 

N0»?««HP 

00002* 

JVDR  *  0 

06002T 

nCOj  *  0 

0  0  0  030 

NSTt 1) 

*  1 

000031 

NST ( 2) 

*  P01 

00203? 

NtD ( 1  ) 

*  <NSET(1)*J)»10*1 

0  00  0  3* 

NED ( 2 ) 

*  »' N S E  T  {  2  )  *  3  )  •  1  0  *  201 

c 

r 

r 

hpt N  SET 

TO  2E»C  To  PRJNT  INITIAL  C^NH  I  T  I  DNS 

000042  RETURN 


00004?  g  mq 

SB 


r  ^  r.  «  m 

f}  A  p  *.  <■, 

oooco? 


0000C2 


m  mi  ^  /  s  |  ■-,  **'  /  jj  /  ■  /  f  /  f  j  /  f  ^  i  *j  j  4 

p  C  "  w  *  *  /  **>»•/  *•  f-  -  J  *• ' 

vDSic*.  ’ 1 . . : f 

dimension  *'i:5n»ii.noe.«ai.Tt(!n,|!i  r c » e*  ?o i  < 3n ; l* ;  ?3 , i? > 

1  ,  r- x-P?  ?)  ,  "*  !  2  s  ,  SmPQT  !  2! ,  C*  (  20.  2  1  .  »S<  20.  2) .  V  f  1 .5,  20 , 2* .  af  T*  f  <  '»  ?> 

2  , VU( 23,?), vpOJT i 20,?) .NSeTf?), g“i : *.2^.2 5, 0“(l 0*20,2) 

X  , fH!NT{20»2) 

4  ,TITL;!2«! 

5  ,TL«10.20.2).  KP!5P(2".?! 

6  ,  iN  I  ( ?>#«i<?0,2) 

7  ,  M  12  ,  ?  !  >  **13  ;  ,  *TA0(  20  5  ,  *?A0<  20  !  ,  T*r  2  )  i**8U  ,35  .  *t  <3  ) 

6  ,  f'va  C  4,  3  ) .  «6D  ( 6 ,  3  > ,  P»4G  T5 » 

E3Ulv*.f N:t  irc'ii,k-)([T*ijj5jyH^i,(Ti(?5n)»*n»(r»MOi, 
10F5T0C  > »  I  TAM79)  #ENMP) ,  <  TA  125*  >  »E-r.i) .  (  T*(50  1  ,FA  > ,  I  T*  (  4«6  1  ,  SUSP?) ) 
2,(TAU6)»$:-n3T>,<TA<56  5,9S)  .  IT  4  <  1  3*  )  .  D* ) ,  (  «D  1 t  *  0 1 )  .  V  )  . 

3  (  T  A  C  6  4  ?  ,3fcTA),  {TA{311)*VC?,  (  T  A  (  4 1 8  >  ,  U )  »  I  T  *  (  2  6  «  )  ,  T  8  )  ,  <  Ta  (  26  5  )  , 

4  T$Tj>«(TA<3e>tf><,  <TAJ266»  .1*0)  .  <74(4.06), EHAf*',  (?Af46)*N?Er)t 

5  (9D(t).6H),  ( t A<5e) , 6*1  V/  )  ,  ( dP  < 401  )  .0*)  .  (Ta(1),T?Yl6) 

6  ,  (  a  U  <  8  0  1 1  >  T  L  )  ,  (  T  A  (  2  7  0  )  j  *njSC>,  I  T  4  (  52  !  ,  D^T  ) 

7,(TA(3*l>.H*nvS).  IT*  S  451  )  ,»w  !  ) , <  x- (13).  !  PL  0  ?>  .  <  T  if  2  0  !  .  T  *S  7) 

6.  (fC(  3  )  .  0<  >  .  ;V*(127.J2),  1^(10,  !  S  Z  >  , ( T  a  <  5*  7  > , °E  1 

9  ,(T4(434),WM03),t7A(4fl6).cDS^P! 


0  0  0  4*  c  ? 
000005 
000007 

030014 

000016 


AX  =  A (  X  ) 

P*J3  ,  . (=NMP(1>  *  FVMp(2)) 

HAN  «  S2»T  (AX/J,  14359! 

C  T  A  *  AX  t  AXHJS 

IF  I  HAT  >OF  ST  30  31.31.3? 


00  0  020 

31 

900*1=)  *  -En*»<15  *  (l,/OF5?Cr 

C  00  025 

FCPH3  •  FaU)  •  11,  *  1,/OFST 

0  00  0  30 

GO  TO  33 

000030 

3? 

9COH3  1  -tXMP<?)  .  (OF6TPC  ♦ 

000034 

rc?M  <  r  a  (  2  )  •  ( 1 ,  ♦"FSTOC  )  / 

030036 

3  1 

9tCO?H  I  HCO-8  *  H*N  /  ET* 

CCOOM 

9E0-'HH0^S«SDSf5D 

0  0  0  0  4  3 

rvAPi*fA(l)/EinOT(l) 

000045 

rvAP2«r*(2)/EHP0T(2) 

00004? 

r  y  ap  «  (FA(ii  .  «=■  a  <  2  )  )/E-0 

000051 

9VAP1  *.iN**3  ( l  )/6-00T  (  1 ) 

000053 

9vAP?*-i\k,5(  ?)  /  E-p7T  (  2) 

0  0  0  3  55 

9  v  A  P  *  3  4 1 P /E -0 

0  0  0  ?  5  6 

!  TEST  .  ”  '  ;  _CT  •  ! P L jt  - 

•  00 „u3 

r F  ( ! T T  j  r  )  i 2 a :  3  24  ,  3  ? 8 

000064 

12  4 

0  C  1*  T  i  N  j  * 

030364 

GO  75  1 

370066 

32* 

*H(  jl  )  *  X 

0  0  C  0  7  r 

A jf  J7  ,  1  )  *  H  A  7 

30007? 

A  k 1 J  2  >  ? )  *  F  VAP 2 

000077 

A  i.  (  .■  1  .  1  )  »  F  V  A  P  1 

00007* 

A  0  <  J  7  '  *  )  *  F  ¥  A  P 

Q  0  0  3  7  * 

A  ^ ( J  L .  3 !  *  9X4  =  2 

000107 

As.,  !  j7,  A)  «  4  X  A  5  1 

030101 

A  w  (  j  l  ,  ?  '  *  H  V  A  » 

0  7  0  1  0  3 

A  J  (  »  7  ,  «  )  *  t  —  *  3  M 

1,1  /  E*?7 

/  E-C¬ 


IO  / 


5  $  q  i  e « 
JOOiOft 
i  .1  o  i  o ? 
0001  ?  1 
o  o  c  i :  ? 

03011* 
0  301  Is* 
OOOli* 
030123 
00Ql2n 

U  0  0 1  2  ? 

33013? 
00013* 
0  0  014ft 
000146 
C00172 

000172 

000206 

0"020ft 

000222 


000222 

00822A 

003226 
000227 
030232 
030233 
00024? 
000253 
000265 
000265 
003267 
0  C  0  277 
000300 
000310 
000312 
000322 
000322 
030324 

000373 

030373 

000376 

000377 

000403 

010403 

030404 

000407 

000410 

000417 

000430 


C 


f 


*j(JZ,  SO)  *  f  CJ*'? 

AU?JZ.19>  *  » 
*U<JZ.?3S  «  *0 
*U«J2.35)  «  TU 
*U ( J2  #  4  3  )  -  fSTG 
JH  ■  J7  *  1 
1  CO'-'T  !  N  J  = 

325  !B7  *  !3Z  ♦  1 


j  L  L !  f  «  MOO( I bZ  -  1,  KPRINT)  .  M E .  1  .AND.  KPRJNT  ,WE, 
IF  ( ILuJT)  GO  to  60 

WRITS  (  6*  330  )  CTITLEU>,  !  •  1*  20 
33C  FORMAT  <1H1/  ( 1M  ,  10A6)) 

wRtTE<6.345)K#U,TG»TSTG,P,P0#RAT»EMACH 
340  FORMATt 3H0X«F6.3,4W  U»F7,1,4H  T«F5,0|4M  TO*F5,0*4H 

l  4H  P0«F7,2,6M  n/F«F7,3.7H  MACW«F5,3> 

WR!TE(6<350)FVAP1,FVAP2,FVAP  »rCOMR 

353  F05M*T(26M  VAPORIZED  FRACTION  0F9.6.4M  FF9.6.7H 

1  #  3X *  4  4 COMB  #  F9 , 6 ) 

WHI T=( ft,360)RVAPl,RVAP2,RVAP  ,RCOMR 
360  FORMAT ( 26m  VAPORIZATION  PATE/JN  0F9.6.4H  FF0.6.7M 


1 


P«F  7 , 2# 


80THF9#6 


0  0  T  M  F  *?  ,  6 


1  ,  3X,4-4:OM3.,  F9.6) 


370  FORMATt70H3OXiniZtB  OROPS  R ( M I L )  MASS(LB)  VtlN/SPC)  U-V/A  TPM 
IP  N jiBER  ) 

60  I  '»  1 

V«NS£T (  I  ) 

HO  350  J«ltN 

vdout t j# i )  *  vntj.n/snsPD 

375  RPD  *  OMtl, J, J }/  EMINTtJ.j) 

37;  FRACT  ■  (EMlNT(JiI)  -  EMtl»J. I >)/EMIN7<J. 1» 

GO  TO  2 

37*  HMM  *  J  *  8 

AJ<JZ«l#MVMJ  ■  «S( J» I >*l,F*3 
mZM  *  J  ♦  28 

A  U  (  JZ  •  1  «  MZ  M  )  •  ABS  (VDOUTtJ.m 
M7N  *  J  +  49 
aj(jz-i.mzn)  *  v(i#j«n 

?  CONTINJs 

IF  ( ILL! T)  GO  TO  380 

*R!TE<6#3?0)J,RS(J,n,EM(l,J,!>,V<l.J,l  ) » VOnuTt J.  »>,Tl(l«J.t). 

i  oN(j,n,  tract*  r?o 

393  FORMAT ( l3X . 1 2i 3PF9 i 4* 0PE13. 5  .  F5 . 0 • F9 , 4 » F8 , 2 i El 1 . 3» 4X , 

1  Ell | 3#  4X -  Ell , 3  > 

390  CONTlNJr 


r 

IF  J  I  Li.  I  T )  GO  TO  61 
uMlTSt  5.400) 

400  FORMAT ( 12H3FUEL  DROPS  1 
61  I  *  2 

M *  >i5 5 T «  I ) 

no  410  j«i#n 

voouTtj.i)  «  vnt j. n/SDSPi) 

405  RPO  «  D4( 1 » J» I) /EMI  NT ( J,  I) 

407  TRACT  ■  (SMINTtJ,!)  •  EM ' f , J, I ) ) /FMI NT { J, !> 


60 


,  -  ••  -  4  »  J 

3  S  ***? 
•**;««« 
033454 
050455 
550465 
J  5  0  4  6  7 
000477 
000477 
000501 

00055* 

000554 

000557 
000560 
000561 
000571 
00057? 
000574 
000576 
0  00603 
000616 
000635 
00065? 
000666 
00)703 
C0C714 
000730 
000746 
000767 

001017 

001024 

001026 

001030 

00103? 

001034 

001036 

001040 

001046 

001056 

001063 

001101 

001124 

001131 

001133 

001146 

001160 

001167 

001177 

001206 

001210 

001217 

001226 

001232 


» ;  *  •  j  *  is 

..  ijM.-f’  «  PM,**  3 

**  { »*  *  J  -  3  * 

. -/*o  *  *bs  <vocj’(w*n) 

-pi  J  *  39 

c  v*l,  j,  I) 

*  CO ^ 7 !N j* 

jr  nun  go  to  4io 

w«!Tr.  <6.391)  J  ,«S  (  J.  P » 1 1 » J*  P  .  « I  1.  J.  P « VOOUT  I  J,  1)  . 

3  TL(  1,  J.  !  > .  0N(  J.  I  > , TRACT,  RPD#  KDtSCIJ.P 
391  C09H4H 13X. !2.3Pf9,4,0PEl3,5  ,  F5 , 0 , F9 , 4 , F8 , 2 , £  1 1 . 3 . 

1  Ell , 3 , 4X , El 1 , 3. 1 15 ) 

410  CO^TINJ: 


no  JO  *  *1  >  2 
no  10  <*1,12 
10  S<*. i ) *3  . 

DO  12  ! *1. 2 
m*njSET(  !  ) 

00  12  J*1.M 

IF ( V  ( 1  J. I ) ,L6,0.  )  60  TO  12 
sii#  n  «s(i»  i  ♦♦dni  j#  p/vii,  j,  n 
s <  2«  P*5(2»  I  ♦♦DNIJ.  I  >/V(l,  j,  I  )*6M(1,  j,  p 
S(  3#  I  )  *5(3,  n*DN(  J,  P/VC  l*  J,  !  >*«S<  J.  ! ) 

5(4*1  )  *  3*4,  n*DN(  J.  I  J,  P*«S<  J,  !  )**2 

5(5# !)*5(5, I > ♦DM ( J. I )/V(l, J, J )*BS< J. I >**3 

5(6,  I)  *5(6,  M*0N(  J,  |  J  J 

S(7,  I  )«S(7,  P*DN(  J.  J  >*OMIl,  J,  P*»S(J*P«*2 

S(6# I >«S<8, I  I ♦DNC J# I >/V(l,J, P-BETAI J, ! >*HS( J,  !  ) 

5 ( 9, I ) *3(9, I l*0N( J# P/V<1,  J,  P*BFTA(  J#  P*RS(  J,  P**i,5 
5(10.  n*S(10#  !)*CD(RH(  J,  !>)*ON(  J#  P*6H(1,  J#  I  ) 

i  /«moit(.(i,  j,  i  >,  n/vu,  j#  n/»s(  j#  p 

J,  2  CONTINJ5 

xm(i>*s<2.i>/su,i> 

XM(2)*S<2#2>/S(1#2> 

XH(3>«(S<2.1»*S(2,2))/(SI1,1)4S<1.?)) 

no  14  i«i#2 

M«MSET(1 ) 

no  13  j«i.n 
vTAB( J)*TU(1. J. 1 ) 

13  XT*0< J>«6H(1, J, 1 > 

CALL  LV1R! VIXHII ),Tm( ! l,XTA8,7TAH,0ynx,N) 

14  RH<1# I)»< .236733  *xh < P /RHO ( T* | p , P >•*, 3333333333 

RHH3»(9-(0(TM(l),l)*(l,*FVAPl)*EH00Ttl  )  ♦H*0<  Tm{  2  )  ,  2  )  *  1 1 , -F 

1  £  <*  OOT  (  2  )  )  /  (  (l,*FVAPl>*E«DOT(l>*(l,-FVAp?>*6HDOT<?>  ) 

R0(1.3>*< ,236733  *XM < 3 ) /PH03 >**, 3333333333 

00  15  ! *1 #  2 

»0(2.  P«(S(5. 1  )/S(l,  1  P**,  3333333333 
P0(3, I ) «SU9T  <S<  5. | J /SI  3, J ) ) 

P0I4, I ) *S ( 5» 1 )/S(4.  1 ) 

R0(5#P*SO9TfS(7,p/S(6,!)> 

RBC6,  P«(S(9,  P/SC8,  !  >  )**2 

15  COMTINJ: 

RBI2.3)«((5(5»1)*S<5.2)>/(S(1.1  )*S( 1.2 )))*♦, 3333333333 
RB(i#3)«S04TMS(5,l)*S(5*2)»/(S(3«l  )*S(3.2M  > 
R0(4,3)*(S(5.1)*S(5,?))/(S(4,P*S(4,2)> 
R8(5.3)*S09T((SI7,1)*S(7.2P/(S(A.1)*S(6.2P) 


o?  >  • 


001240 

0012'- 
0  0  S  246 
0  j  i  2  4  7 
001251 
00175? 
noi?54 
00125* 
U0126P 
001270 
001274 
001276 


001322 

001324 

00132ft 

001327 

001342 

00134ft 

001347 

001351 

001353 

001353 

001355 

001357 

001360 

001362 

001363 

001366 

001370 

001373 

0  013  7  3 
001375 
0  01400 


001400 

051430 


0  0  1431 
00145? 


00115? 

001460 

001460 

001465 

001467 

001471 

62 


9&<6,3)*<(3(9,l!*S(9,2))/(St8,l)*S(9,2>)>»*2 

r 

RANE--RAN/ETA 
XL<H«9VAPi*RANE 
XL< 2>*9VAP?*RANb 
XLI3)*9VAP  *9A\E 
no  17  Ill,? 

M*VSET( S ) 

DC  17  <*1,6 

CALL  LM3«l'MRrt<K,  J  )  ,DVfl(rf,  |  > ,  RS ,  VOOUT  ,  rs  Ynx ,  N ) 

17  COVTlNJs 

no  19  <*1.5 

nv?<<,3)*UBS(nvH(K,i)>»(i,-rvAPi>«EMDr'T(iu 

1  ABS!DVfl(«,2))»(i,.rvAP2)«EMDCT(2))/ 

2  Ui  s  »FV  API  )*EMDOT(  1)*(  1 ,  •FVAP?)«EM[pT{  2)  ) 

18  COMTJNJ; 

no  19  1*1.3 
no  19  <*1-5 

19  9EP(<.  n*H3(X#  n*PEC 

nC*e  375°PA^E*fc7Ar'6X 

ORaGU*  *  StiO.D*  DC 

09*0(2)  =  S ( 10 1 2  )  *  DC 

l—fcGl 3  )  «  (S(lOil)*  S ( 10 , ?  )  )  *  DC 

c 

GO  TO  4 

20  AU(J2-l»82)  *  960(6,1) 

AUUZ-1.8j)  «  9£n<6,2) 

AUUZ-1,84)  *  RED ( 6,  3 ) 

AUtJZ-1.85)  «  9tC0*B 
AU ( JZ” t «  86 )  *  XL  <  3 ) 

All  f  vZ”  l »  97 }  *  ABS(0«B(6,1) ) 

A U  <  JZ" l , 88 )  «  ABS ( UV8 ( 6 , 2 )  > 

A U ( JZ“1 ,89)  *  A8S(DV0(6,3) > 

4  COvTINj: 
r 

If  ( I LL I T )  GO  TO  62 
WRITE  (6,90) 

9G  r09MAT(lH0,///8X,6MPADIUS,l5*.8HOX!D!ZER,;2X.4wrUEl ,26X,5hT0TaL  ,/ 

1  24Xt iH9,8x,2Nnv,8X,3HRE0,»X,lM»,8X.2KDV,8X,3HPED,8X ,1MM,BX, 

2  2MDV #  9X , 3WR6D, / /  ) 

WRITE  (5,91)  ( (R8(K, I ) ,DV8(K, I ) ,RE0(K,  I  > .  I  *  1,3),<  »1,6) 

91  rORMAT(20H  MASS  mean  D80P  , 3 C JPF 1 0 , 4 , OPf 1 0 . 4 , r 1 0 , 1 ) , // , 

1  20H  VOL,  mean  DPOP  , 3(3PF10 ,4, OPPIO , 4,H  0  , 1 ) ,  /  /  , 

2  20H  VOL,  OJAM,  MEAN  , 3 ( 3PF 1 0 , 4 , 0 PF 1 0 , 4 , F 1 0 , 1 ) , / / , 

3  20H  VOL,  SUPFACF  MEAN  ,  3  < 3=>Fl 0 , 4 , 0PF1 0 . 4 ,  FI t ,  1 ) , // , 

i  2CH  MASS  LOSS  MEAN  , 3 (3PF10 . 4 . 0PF10 , 4 , FI 0 , 1 ) . // « 

5  20H  STABILITY  MEAN  , 3 < 3PF1 0 , 4 , 0PF10 , 4 , FI  0  . 1 > , / /  ) 

WRITE  (6,92)  <XL(i>,I«i,3),9LC0M8,<DPAG(!),i»l,3> 

9?  roQMAT(  ///,31<.8H0XIDIZEP,4X,4HFUEL,*);,5NT0TALh  4X.8HC0MBUST  ,  ,/ 

1  ,23-t  BJ.9NING  PATE  P  AR  AMETEP.  t*  ■  '-r  1 0 , 6 ,  / 

2  ,15M  DRAG  PARAMETER,  ISX.SFICS,  2,/} 

WRI TEC  6, 56)DX 

55  -OPMAT(19HO!NTEGHaTION  STEP*F7.4) 

C 

62  nA»RAN/9H03/SDSPn/AX/,592 
DAOX*D*»EMOOT(1 )*RVAP1 
OAFU«UA*EMOOT(2)*RVAP2 
nATOT*OA*(=MDOT(l)*EMDOT(; : )*RVAP 


0  0 1 4  7  < 
05147? 
001903 
061906 
001914 
061321 
001523 
001549 
061344 
051947 


RACU*D**«S4D0TC1UFmD0T(2> 

*S}0X*»a**R9a»1«6NR0T<1>/AX  *1 ,4349/1, 3^01 

S.RrU*R*'4*HVA42*EHD0T(2)/4X  *1  ,43P9/l,u«Ol 

HSr?0«R*  *J?*DOT  <  1  )*gHPftT  (2  > )  /  AX* 1 , 4389/1 , 0 60 1 

«RC0«'.*X*«:0^H*<!;f460T(i}*E‘‘0i'Tf  2»  )  /  A  X*1 ,  <399/1 , 0401 
;r  (iljt)  go  to  63 

w« T  T  =  <  S ^  93 )  nAOX,bAru,OATOTsnArOjb'}OXfH^fU«BRTr)t^»rC 
93  r0RM*T<///.3OXt4FlO,2,/> 

63  RETURN 
fcNO 


63 


C230SS 

OOOGO* 

050003 


000003 


00000^ 


050011 

1*00513 

00001ft 

080020 

050022 

000024 

000025 

000027 

00003ft 


00003ft 


030041 

000055 

000053 


000054 

limin') 
*,  v  *,  o  v  . 

050064 

Q30071 

00007? 

0000/7 

000101 

000101 


S jPHDUT I Nfc  PRGP(NSP) 

COHMQM  /E/re.MH/e/8C/C/TA 

DIMENSION  PC(«li),  NH  <  20 ) ,  BD ( 10  #  20.  1?). 14(  1000  ) 

ECU1VA.ENCE  ( T* { 30 ) »  P),  { T  ^  s  25  ? : .  rat;,  (  TA<2ft5),  TSTG) 

1.  ( T  4 { ?64  ) ,  TG>,  <T4<2ft6)#P0).  (TA(2ft7>,  V!S>.  <TA<2ft8),  F*G> 

?,  CP3),  t  TA  i  39i } ,  Rm  ovs>j  <lM4a3>,  WTML» 

3,  (  T  A  C  4  34  7 »  RHOG) ,  <TA{ftB5) ,  eMAC"),  (  T  A  (  «0ft )  4  SOSPD 1 » (  T*  I  ?*« )  #  GAMl  ) 

4,  (TA( ?ft0) ,G4M2), (T4(26D .R4M35 » ( f  A  (  2ft?  ?  » GAMA ) » <  T  A ! 2A3 ) , G A^S ) 

5» ( TA ( 25) *  H ) »  (TA(27>,  R ) ,  f  TA { 47B ) ,  m 

0 

C 

C  *****  33T4IN  STAGNATION  TEMP, TSTG,  AND  T«E  DFRIVAMVE  Of'*  TFM® 

r  with  respect  to  *vf,  trod 
r 

CALL  TGiNT  (BATjP.TSTfj.TGOD) 

r. 

C  *****  OBTAIN  MOLECULAR  WFJRMT  OF  PRODUCTS  AS  FUNCTION  O/F 

CALL  HT4 1  NT  <  R ftT ,  P ,  WTHL ,  WT HI.  DT  > 

C 

RG*H/WT'1L 

3A4  =  F  GAM  <  °*7  ) 
r,AMlsGA<1-l, 

GAM2«G411/2, 

GAMJ*GA  42/ 3AM 
RAMfttGANl/SAM 
GAM5  *  GAM/GAMI 
U2»U»U 

c  *****  CALCULATE  static  ten®,  ►  ROM  TOTAL  TFMP,  and  RAMMA 

TG»TST5*<1,-5AM3*u?/RG/T^T(./G) 

c 

c 

c 

r  ;**•*  CALCULATE  GAS  DENS  I TY 

rhog»p/rg/tg 

t? 

r 

c,  *****  CALC  .".ATE  VlSCUSITV,  TMERMAL  CnsjOUC T  !  w  \  r  v ,  *vd  S°EC1F!C  wfc*r 
r  OF  GAS  AS  HJNCTiON  OF  P,  TO,  O'T 

CALL  VtSG.'S  ”SAT,P#TG.  VIS.FKG.CPG.NSP) 

C 

PHnys-i^-ic 3/v is  2, 

RHOTF* , J75*RW0G 
r. 

r.  V.CILATE  T0?AL  PRESSURE 
»0  *  P*(TSTG/TG)**GAM5 

)C  S  Kj  G 

SPSPO  «  5UTT (RR  *  GAM  *  G) 

FMACH« J/SUSPP 
IF ( EMAC4  ,3T,  1,0)  STOP 
EM. Ci2  »  t  .  A  BH»  •*  2 
C 

?E  T JRN 
ENT 
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SUBROUTINE  REED 

30360*  DIMENSION  *e<80>#  NM!?!)),  BD(10,20,12>,TA<in00>,T!TLE(2<*> 

00000?  COMMON  /NLN/  LINES*  LNS 

000002  COMMON  /E/  F  C ,  NM  /8/«D  /C/TA 

000002  COMMON  /  ADDPPM  /  PARI,  P  AR2  •  PAR3 »  PARA,  PARS  *  P AP6 

00000?  COMMON  /PRNT/  mPRJNT 

000002  COMMON  /PR3S/  PRPSYS<4)  /SPEC/  NSp,  SPNAM|(2,  15)  /O F/  Dlll))i  nr 

1/QTL/  NTL<?)  /DTfLR/  NTFLM(2)  /8T/  NT  /PRSSR/  PR$$(3)  /CMPTAR/  VSP 
2(15,  15,  3)  /NTM/  RMTaR(15,  3)  /TQ | /TCI C 15#  3)  /TBFGAM/  FGAMT (15) 
3/TA8T/  TTA3(I5)  /V*SGS/  VSP(15.  15),  CPTAfHlS,  15)  /MR/  uM18) 
4/DIFFuS/  E»S)(15;  2).  S!GNa(15,  2),  EPS0F<2>,  SIGOr(?)  /TL*/  TL i 15 
5,  ?)  /»**/  P*(15,  2)  /RHOT/  RmOl  (15.  ?>  / CPr/  CP(15,  2)  / 1  *  M  7  / 

6L AN <  15 ,  2)  /T RLMX/  TFILMU8,  2)  /VISCV*/  WIS(15,  2)  /*VAPV/  KVAP 
7(15,  2)  /C»VPY/  CPVAP<15,  ?>  /LINKA/  LOOSE 
C  00002  EQUIVALENCE  (TA(1),  T  J  TLE  ( 1  )  ) ,  <TA(25),RR).  (Ta(26),  TWOPI) 

1,  (TA<27),  G),  ( T  A ( 26 ) ,  FTP!),  (TA<29),  ThRR),  <FC(t).  P 

2,  (  TC  <  6).  UMAX)  ,  <Fe<R),E*IN>,  ( PC ( 9 ) «  EMAX ) ,  (FC( 14 ) , T^TOP* 

3,  C  NM  (  7 ) ,  NSW)  ,  <  NM  (  9 ) ,  “P),  (NM(H),  NO)  ,  (NM(lO).NV) 

4,  (TA(30).  P)  ,  ( T  A ( 31 ) ,  RATMN),  ( T A ( 32 ) , RATmX  ) ,  < T A ( 33) , AX* I N ) 

5,  ( T A ( 34 ) ,  AXMAXJ,  (TA(35),  OFSTOC)  . C T A ( 36 ) .  V)  , ( T A ( 36 ) , ENflOT ) 

6,  C  T  A  <  40  >  .  TLL  )  ,  04(42),  RM  ),  (  T  A  (  44  )  ,  SJG).  (TA(46).  N*fcT) 

7,  (TA(48>,  mTMOL ) ,  <FC( 11) ,  W>  ,(TA(478),  U)  , < T A ( 412 ) , 0 JR ) 

8,  (  NM  1 5 ) ,  NPLOT  ) 

00C  30?  u  l  MENS  I  ON  V  <  2  > ,  EHOQT  t  ? )  ,  TLL  <  2  )  ,  RM2,.  mh<?>,  NSET(2) 

1,  4TNUi.<2)  ,  w  ( 3 ) 

C00002  RR*  18540, 

000004  TwOPl  s  6,2881853 

000005  G  «  386,4 

0  00007  FT 0 !  a  4,1887902 

000010  THRU  »  ,33338333 

000012  G JR  «  8.78E-10 

000013  READ  (5,  1)  (Tl)LE<I>»  I  •  1,  24) 

000025  1  FORMAT  (1046) 

000C25  WRITE  (  5,  700  )  (TJTlEP),  f  •  1,  24) 

000037  700  TORM4T  (1H1.10A6/  (1X1«AM) 

000037  RE  4  3  <  5 , 5  )  T,  HWAX,  Emin,  E«ax,  TSTOP,  u 

000057  5  f 0 R 4 A T  (7ElO,b) 

000C57  READ  (5.8)  NSW »  MP,  NO,  NPLO.,  KPRJNT 

000078  8  rOR4*T(6!10) 

000075  READ  (5.9)  P,  RATmn,  RATmX,  AXMJN,  AXM4X,  ofstoc 

000115  RE  4  0  (5,9)  (EmDOTU),  V(1),Tll(!),  R*  (  I  ) .  S  I  G<  H  ,  N5ET  <  I  ) ,  I  «1 ,  ?  ) 

0  001 46  9  F0R<<4T<5E10 ,6,  110) 

000146  RE  *0  (  5,  6000  )  PARI,  PAH2,  P4R3  PAR4,  PAR5,  P4RA 

000166  6000  FORMAT  ( 6F 9 , 0  ) 

000166  w  ( 1 )  *  1,0 

000173  *<2>  *  1,0 

030171  *<J)  *  1,0 

000172  N M 1 3  )  *  1 

000173  NM  1 2 )  •  1 

0  001  74  NM  14)*1 

000175  WRITE  (6,  801)  PRPSrS 

0  00702  602  FORMAT  (16H0!N®UT  F QH  OX  a  2A6,  9h,  FjFl  *  2*6/) 

30020?  WHITE  (5,  70)  NH ,  NS°»  NTl ,  NTF^M,  NT 

000220  70  F0PM4T  (  74H  NU“BER  OF  O/r  R4T|*s  .  \?/  21h  nu- 

1REP  OF  SpE:iES  •  I?/  38M  numNER  OF  LI3>'IP  TpMpfp*TiRES  (OX)  •  12/ 

2 4 0 "A  NU«36M  Or  i. I OU 1 0  TEMPERATURES  (f  JED  »  !?/  36*  mumocp  of  F  1 1.“ 
3TE«PEP*  TUNES  (OX)  «  52/  38h  NJ-8FR  Or  fjlm  rt-Pf-D » TyoE*  (FUEL )  •  I 
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42/  36H  NUM9ER  Of  GAS  PHASE  TEMPERATURES  *  I?/) 

C90220  LINES  »  14 

090221  NTll  *  >4TL  <  1  > 

000223  NTL2  «  XTL C  2  > 

000224  NTFL^l  «  NTFlM<1) 

09022A  NTFLM2  ■  N TF u M < 2 ) 

090227  WRITE  (5,  701) 

000233  701  FORMAT  <IH  O/F  RATIOS/) 

000233  WRITE  (6*  702 )  <R(!),  I  •  1,  NR ) 

000246  702  FORMAT  < 12X6E20 , ft ) 

090246  IF  (NR  , GT ,  12)  LNS  «  5 

000252  IF  (NR  ,LE,  12)  LN$  «  4 

090255  IF  (NR  ,L£,  6)  LNS  «  3 

030260  LINES  •  LINES  ♦  LNS 

090262  00  15  1  NO* 1 #  3 

000263  WRITE  (Si  70S)  INQi  PRSS  UNO) 

000273  703  FORMAT  (4M0PC(lli  4H)  •  Elft.ft,  5H  PSIA/) 

000273  LINES  *  LINES  ♦  J 

Miiii'ri  IF  (i.iv  =  S  ,G2,  5C)  WRITE  (A,  704  5 

090303  704  FORMAT  (1H1) 

090303  IF  (LINES  ,Cfc,  50)  LINES  •  0 

003307  WRITE  ( 6 i  705) 

090313  705  FORMAT  ( 23 H  SPECIES  MOLE  FRACTIONS/) 

000313  LINES  *  LINES  ♦  1 

000315  IF  (LIN=F  ,GP,  50)  WRITE  (A,  704) 

000323  IF  (LINES  ,G6,  50)  LINES  ■  0 

090327  00  12  I  *  1  *  NSP 

000331  WRITE  <  5.  7020  )  SPNAME < 1 •  I),  SPNAHE(2i  I),  <YSP(J,  I.  IND^i  J  *  1 

1.  MR ) 

0 0 0 35 A  7020  FORMAT  (1H02A6,  FlV.ft,  5F20,8/  ( 12*6F20 , ft ) ) 

09035A  IF  (NR  , GT ,  12)  LNS  «  4 

00036?  IF  (NR  ,LE,  12)  LNS  "  3 

030365  IF  (NR  .LE,  6)  LNS  «  2 

39037C  LINES  «  LINES  ♦  LNS 

000372  IF  (LINES  ,GE,  50)  WRITE  (A,  704! 

0  00  4  0 11  15  IF  (LIMES  i  GE  •  50)  LINES  «  0 

000407  WRITE  (J.  706) 

090412  706  FORMAT  (18R0AVE,  ^'S  MOL.  WT,/> 

330412  WRITS  (J,  702)  <wmT*8(J,  inQ),  I  »  1.  NR) 

090*27  CAlL  P A  jES ( NR  ) 

000431  WRITE  (6.  707) 

000435  70  7  FORMAT  ( 32R0EQU  I  LI  BR I U*  fi AS  TE*P,  (RANKINE)/) 

090435  WHITE  (6#  7^2)  (Tci(I,  IND).  I  •  li  NR) 

000457  15  CALL  P A  j£R ( NR  ) 

00045*  WRITE  (5.  70ft ) 

000  46i  704  FORMAT  (6H0GAMMA/) 

000462  WRITE  (6.  7  q  2 )  (FGamT(J),  j  •  1,  NR) 

000475  CA„„  PA2IR(f;t) 

0  0  C  4  7  7  WRITS  (6,  7  0  R ) 

090503  739  FORMAT  ( 26RQGAS  phase  TE*P.  ( R A Nk I Nf  )  /  ) 

000503  WRITE  (6,  7Q2)  (TTAy(j),  j  ■  i,  *T) 

00051*  CALL  P A  jER ( NT  ) 

000520  WRITE  ,‘i,  710) 

030524  719  FORMAT  (5H0GAS  PHASE  V1SC.  (  (L9S  «ASS  /  IN  -  SFC)  X  18  •*  *6)/) 

090524  LINES  «  LINES  *  3 

00052A  IF  (lINsS  ,GF,  50)  WR  J  TE  (A,  704  ) 

00053*  IF  (tlNES  ,GF,  50)  LINES  *  0 

000540  DO  17  !»1,NS® 

35 


000542 

WHITE  (5,  ’02)  <VS*><J,  I  >  t  J  *  1.  *n 

000554 

17 

0 ALL  Pj*<NT) 

000543 

write  (6*  ?it) 

00056* 

711 

FORMAT  ( 34  HC  GAS  PRASE  CP  (f»Ttl  /  M-’JLg  •  DeU  B)/> 

000564 

LlVfcS  *  LIMES  *  3 

000570 

IF  (LIMES  .  IjF  •  50)  WB  i  T6  (a,  7*04) 

0  0  0  5  7  6 

I F  (LlMfS  ,i,F,  50)  LINES  *  0 

000*0? 

no  ip  (»i.'i>p 

0*0*04 

writs  a,  7r ? )  (cp?aB(j.  i ; ,  j  ■  i#  NT) 

000420 

1^ 

CALL  P3*(NT) 

000*25 

WRITS  (6.  712 ) 

000*30 

712 

FORMAT  ( 1 7 RO  SPEC T £ S  MOL.  WT./) 

000*30 

WRITS  <  6*  702  )  (WM(  ! ) ,  I  i  1,  MS»> 

000*43 

CALL  PAjEH(NSP) 

000*45 

WRITE  (4.  713) 

000*51 

713 

FORMAT  (35-tOSPECIES  IMTERACTIO’"  F\=WGY  (HER  *)/> 

000*51 

WRITE  ( 6<  702  )  ( EPS  I  (  |  #  1).  I  *  1,  MSP) 

000664 

CALL  P A  jES ( NSP ) 

000*66 

WRITE  ( 6 >  714) 

00047? 

714 

FORMAT  (33wn»K0P,  INTERACTIQM  F  M  F  W  G  Y  (DEC,  K  )  /  ) 

00067? 

WRIT?  (6,  702)  EPSCF 

000700 

CALL  P * 3EP ( 2 ) 

000702 

WRITE  i 6  <  715) 

00070<l- 

715 

FORMAT  (  30  HO  SPEC  I  £S  *01,  u  j  *  ,  (  AMijRTROMS  )/ ) 

00070* 

«R1T5  ( 6 1  70?)  (SIGMaU,  1),  !  t  l.  MSP) 

00C721 

CALL  P  *  3ER ! MSP ) 

000723 

WRITE  (6.  716) 

000727 

716 

FORMAT  <  2tHCaR0p  ,  “OL.  OIA,  (  AMOSTdcmS  >  /  ) 

C00’2T 

WRITS  (6,  702)  SIGCF 

000735 

CALL  P » 3EH ( 2 • 

000737 

WRITE  <*.  717) 

0  0  0  74  * 

717 

FORMAT  (15W0PRCP.  mql,  WT,/> 

000743 

WRITS  <5.  702)  wTmOl 

000751 

CALL  PA  3FR  <  2 ) 

000753 

WRITE  <5.  714) 

000757 

71* 

FORMAT  (26-fOOX,  L  I  0  .  TEMP,  (  R  A)'K  IMF)/) 

000757 

WRITS  (6.  702)  (Tl<J.  1),  J  *  1,  NTL1) 

00077? 

CALL  Pa3ER(NTL1) 

0  C  0  774 

WRITS  ( 5 »  720  ) 

001000 

720 

roRMAT  (25W0CX.  Va®OUR  PRE«S,  ' ® S I  a ) / ) 

0*1000 

wRHS  (6.  70P)  (PV(J.  1).  J  *  1,  mtlI) 

001013 

CALL  PA3EN(NTL1) 

001015 

wR  t  T  S  (  5.  722  ) 

001021 

72? 

r  Q  RM  A  T  {.liiwOOX,  L  1 3 ,  USMSITY  (|RS  “ASS  /  CU  IM^N)/! 

001021 

wRITS  (6.  70?)  <»hPl<J»  1).  J  *  1.  MTL1) 

00103* 

CAUL  P  4  jE  M ( N  1  L 1  ) 

00103* 

WRITS  (5.  724) 

00104? 

724 

tqrmaT  (  36  -AO  OX  .  L  I  3  .  C°  TRTJ  {  uH  “*SS  -  !ER  R  )  /  ) 

00104? 

WRITS  (5.  7C2>  (CP'J#  1).  J  *  1.  N T l 1 ) 

001055 

CALL  PA3ER  (MTL15 

001057 

WRITS  (S.  726) 

'*1067 

726 

roRMAT  (33R0OX,  MEAT  Or  VAP,  ;wru  /  Lg  MASS)/) 

001063 

writs  (6.  70?)  (LA“(J,  1),  J  *  1#  RTLl) 

00107* 

CALL  Pa jER( nT|_1  ) 

001100 

WRITS  <5.  719) 

00110* 

7i? 

FORMAT  (  26  -10 r  UtL  LIO,  Tfcwp,  (3*\<[vf)/) 

0 0 11 c  * 

wRITf  (6#  70?)  (Tl(J,  ?),  .<  =  1,  mtlp) 

09111' 

C a _L  »a3ER(MTl2) 

tmiji 

mi** 

081133 

001140 

00114? 

001146 

001146 

001161 

30116-5 

001167 

00116? 

30120? 

001204 

00U10 

001210 

001221 

001225 

001211 

001231 

001244 

001246 

00125? 

00125? 

001265 

001267 

001271 

301271 

3G1S06 

001310 

00131* 

001314 

001327 

001331 

001335 

001335 

00135C 

00135? 

001356 

001356 

001371 

00137J 

001377 

001377 
o  o  i  4 1  ? 
00141* 
0  01*20 
001*20 
j:i*33 
001*33 
001**2 
001*50 
001*50 
001*50 
001*6? 
08146? 
091505 


WRITg  (6i  721) 

721  FORMAT  (2CH0PUSI  V*»OUS  PRFSS,  t  PS  f  A » / ^ 

WRITE  (5.  702)  (PV<J»  2).  J  •  t»  M?L?> 

CALL  P A  3ER  <  NYL2 ) 

WRIT?  (6.  723 > 

723  Ff«4»T  (J9R0FU6L  LlQ,  DENSITY  (LBS  MASS  /  CU  INCH)/) 

Mfc «  i .  C,  70*)  (WHfUJ.  2).  J  *  %>  NT12) 

CALL  P43SH(WTL2) 

WRITE  <8»  729) 

723  FORMAT  (37HCPUEL  L ! 0 •  CP  t n T U  /  LB  MASS  »  l>E«  Ri?) 

WRITS  <S,  702)  <CPiJ,  ?>.  J  «  t,  N7L2) 

CALL  M3SRINYL2) 

WRITE  <5.  72?) 

72?  FORMAT  ( 34R0FUEL  HEAT  OF  VAP,  (BTU  /  L?  »**SS ) / » 

WRITE  (  8.  702)  (LAMUi  2).  J  *  1.  *?L?> 

CALL  PA3ERIMTL2) 

WRITE  ( 6 »  72B) 

72*  FORMAT  (27R0DX,  VAPOUR  TEMP.  (RANKJNE)/) 

WRITE  (5*  702)  <  TF 1 LM  ( J»  1).  J  •  1.  NTFLHi ) 

CALL  PS3ER(N?FLM1) 

WRIT?  (8,  730) 

730  FORMAT  (52WOOX.  VALOUR  VISC,  <  f  L « S  MAS”?  /  IN  -  5EO  X  10  **  *6)/) 
WRIT?  <8«  70 P)  ( V V  I S ( J i  1).  J  •  1*  ntflmij 

CALL  P  A 3£R  (  %'  f F L M 1 ) 

WRITS  (8,  731) 

732  FORMAT  (62W0OX,  VAPOUR  CONDUCT  I V I TY  MRT'r  /  IN  *  SEC  -  DEG  R )  7  10 
1  **  ♦ ft > / i 

WRITS  (  8.  702)  ( K  V  *  P ( J  *  1).  J  *  It  NTFLM15 
CALL  PA3FRINTFLM1) 

WRITS  (8.  734) 

734  FORMAT  <38-»0nx,  VAPOUR  CP  (BTU  /  L«  MASS  -  DCG  *)/> 

WHITS  <  8.  702)  ( CPV  AP  ( J  i  Hi  J  ■  1.  NTFL*1> 

CALL  PA3ERCNTFLM1) 

WRtT?  (8.  729) 

729  FORMAT  (28R0FUFL  V*POU»  T-gwp,  ( R AN*  JNE 1  /  ) 

WRIT?  (  8,  702)  ’TFILM(J,  2V,  J  *  i,  NTF  LM2 ) 

CALL  F436R(NTFlm?) 

WRITE  (8.  731) 

731  FORMAT  (53-<0S‘UEL  VAPOUR  VISC.  («LBS  M/SS  /IN  -  SfeC)  X  10  •*  *6)/) 
WRITE  (Si  70  2)  <  v  y  J  S  t  J ,  2),  J  «  1*  NTfiH?) 

CALL  P43ER(NTFlM2) 

WRIT?  t  > .  733) 

73  5  FORMAT  (63R0FUEL  VAPOUR  CONDUCTIVITY  ( < H Tu  /  IN  -  *EC  -  DEC  R)  X  1 
10  **  ♦A)/) 

WRITS  (Si  702)  i*V*PU*  2),  J  *  1.  NTFl.M?  ‘ 

CALL  PA3ER(NTFuM?l 
WRITE  (8.  735) 

7JS  FORMAT  (39R0FUEL  VAPOU»  CP  (BTU  /  iB  MASS  -  0£S  »)/) 

WRITS  (  8.  702  )  '■  r.?v*P  (  jt  2).  J  •  1.  N  T  F  L  “  i' 1 
CALL  P»16R(NTFlm?) 
i F  (wOOSE  - €0 ,  0)  WRITF  (6,  736) 

IF  (u003f  ,fc3.  II  MRJTE  (6,  737) 

736  FORMAT  (24W0WO  Fuel  TwERhal  OfcCOMP,) 

73?  FORMAT  (2lwfiFU6L  THERMAL  DFCCHP.) 

16  WRITS  (8.  ?0)  (TiTLEUl*  I  *  1.  24) 

20  FORMAT  (44R11SN7002  S  TE«  0  Y.  ST  A  T  E  SP9*Y  COMBUSTION  »00£L/<  l4" .  10*6  D 
WRI T?( S,*0> (5MROT( I ),RM( I) ,S55( I ).*TMSL U  > .  1*1.21 

<0  ruRMAKiHO 


pFTf 


Q  « 150  *> 

P0151S 

051115 

311514 


1  8*  ?Cm  M A 5S  FLOM  RH  ■,» 

0  F"  1 5  5  2 «  El?, 3,  ?F7, 

**"  F  FIS. 2.  FI?  3  ?F7 

WRITE  (5,  30  >  MR  * 

50  FORMAT  (19M0MMAX  /  HMJN  «  ?  **.  ,2> 

»?rtm 
fc  vo 


sta 

?F7,2 

?r7.2 


*OlE  XT 
/ 
> 


010902 

09000? 

090002 


09000? 


ooooc? 

000003 
090004 
000005 
09000? 
0900  11 
090014 
090015 
090  020 
0  90024 
090024 
09003? 
9000*3 
0  00  044 
090059 
09005? 
090069 
030066 
00007* 
030101 
0  90  1  0  m 
000113 
000117 
000123 
090127 
000133 
090134 
990137 
0001*6 
9901*6 


SUBROUTINE  RSgT 

DIMENSION  FC ( SO ) ,  NM<20),  RD( 10,  20#  12)*  T A 1 1000 ) ,F * ( 2 ) . fc^MP ( ? . 
COMMON  /=/  FC#  NM  /»/  RO  /C/Ti 

E3U1VA.ENCE  (NM{141#  JPZ)  ,<NM(1?>,  JZ),  (TA(50),  Fa) 

1, (14(472).  ENM® ) i ( Y  A ( 481 1 #  RW \) , ( T A ( 258 ) . EHD > . ( T A ( 36 ) »  FMD4T  ) 

2, (TA(52).  ONT  >,  (TA(?4),  wCON),  (TA(4?>.  »M),  <TA(44),SIQ> 

3# ( TA( 46 ) #  VS6T),  (Ta(40),TTL>.  ( T  a  « 3«  > .  VV)  ,(«D(B01).  TL) 

4  # (40(1631)*  V)  .  ( 8D ( 1 )  i  EM)  ,  (Ta<54),  EM  I  NT ) , ( TA i 96 ) . «S ) 

5# ( TA( 136) »  DN)»  ( T A ( 1 76 ) »  K  K )  ,  ( T A (?14 ) ,NU ) , ( T A ( 25 ) ,  5) 

6,  <TA(?57).  RAT),  (TA<256),  !)NT  T  ) ,  (  T  A  (  4#  > .  WTMOL  > 

7,  (U(U:,SATMN>#(TA(32)#RATMX) 

njMEvSIOv  fcMDOT (21#  0*1(2).  WCOM?>,  BM<2i.  SIG(2),  n5ET(2> 

1,  TTl<2>.  V  V  (  2  )  #  EMS ( 20  > ,  DN$(20).  RSS(20>.  TL  ( 1 0 . 20  ,  ?  >  ♦  *>“  <  *0 , 2 > 

2,  <K  !  20  »  2  >  i  BS ( 20 , 2  ) ,  DN(29,2),  E“ I  NT ( 20 . ? ) . EM ( 10 # 29 . ? ) 

3,  mo. 20.21  .xtmou?) 

HZ  *  0 

J?  *  1 
PHI  «  0,0 

EMO  *  EMDOT(l)  ♦  F Mf)QT ( ? ) 

DO  100  I*  1.2 
F A (  1 1  «  EMOOT ( I » 

EVMP( I )»  0,0 
D VT (  I  )  «  0,0 

W-3N ( 1  )  •  6,2831853*  uTttOLM)/  R 
N  ■  NRET(I) 

CALL  L33NRM(RM(I),Sir. (I).n.RSS) 

CALL  MASSET(FMS,RhO(TTL( !  1 , J > . HSS . EMDOT ( I  1 . N , PNS > 

NSETU)  •  M«* 

N  «  vS5T(  I  ) 

nc  loo  j* l , n 

T„(l, J, !  )  *  T T L <  1  > 

V( 1. J.  I  1  *  V V ( J  ) 

fc m ( 1 »  J »  I  )  *  E*S< J) 

F  m  |  N  f  (J.u»  EMS(J) 

HS( j, 1 )  «  BSS( J) 

l)V(J,I>  1  DNS(J) 

0 v T (  1  )  *  QNT(  I  )  *  DNS  ( J  1 

K  < ( J, I  )  *  1 

w  * ( J  •  I  )  »  0.0 

100  CDvTHufc 

Osf  T  I  Pint  ( 1 )  *  ONT  (  2  ) 

HAT  «  EMDOT ( 1 ) /FmHDT ( ?  ) 
raTianinH’atmx.ahaxk  RaTmai.batu 
U- TUHM 

fev3 


t 


69 


rUMCT 1  OH  MT(RS) 

C  CALCULATES  DISTANCE  OF  DISSOCIATION  FLAME  AS  FUNCTION  OF  RADIUS 
000003  DIMENSION  RRB  ( 5 ) » RRT < 5 ) 

000003  OATm  RRS,  SR?/  0,,  100,,  200.,  300.,  500.,  13,,  1?,,  1*..  18./ 

C  CONVERT  «S C K »  TO  MICRONS 

000003  RST  *  RS 

0-500  94  RST  «R?T  *  2,54E*4 

00000*  CALL  LN3RIN  <RST,RT,RRS,RRT,DRTDRS,5> 

00001?  RT  •  RT  /  2,54£*4 

000014  RETURN 

000014  end 


■swanupwiwi 


ouhrout  r t  shift 


00090? 

OlHrMSHN  e-C  t  80  )  #  HO  <  2*00  )  ,  NST  (  2  ) .  vEC  (  ?  ) 

(loono? 

/>-'  fc,  5m  /a/  hi; 

0  (J  0  0  0  < 

C.MlVA.i  c=  (NH(5)  ,MPT»0  »  <  *«*4 1  4  1  ,H(,AH),  (  MH  (  3  ;  .  u  A  L  P  >  » 

1 

(MM(l).lH),  (  r  C  <  3  > » H )  1  tfr<  7} ,  -<7D?  )  .  (NMUOl.'JV).  (rCMI,  O), 

2 

( rc  0 1 ) . T  0 ,  ( NH ( 2 ) j  !Nf)D),  (70110),  fl),  <»M(8)»*e«U), 

3 

(  f  C  f  6  ) »  HHAX),  (  NM  (  6  ) ,  MPT  5  ) ,  <rrC5),M“lH),  fr-,i>),  F“*  X ) » 

4 

(t C< 2) #H<T )  , <8M< 17 > ,Nftf ) , f  M*( 195 .  ngft) 

c 

c 

>  * !  F  T  j>l'*TtS  THE  v*^M*8LE5  50  That  o*oppp  POS  I  T  !  O*  I  V(i 

r 

rrCUPS  OUHIMG  T Hf  PROPER  C*CLE 

t  '  r  n  o  ? 

r 

(f  (  M9TN)  15,  10,  1% 

000005 

10 

h?T\i  t  HP TS 

0<*C305 

15 

H*,*-  <  •  '1  0  AH 

0  0  C  0  0  * 

[F  (  HiilH  )  87,  20,  60 

0  ?!  0  0 1C 

20 

[M  H*.P  -  2)  25,  48,  4ft 

030013 

25 

I  “  «  t  4  *  1 

010015 

mPTn  ■  hPT*  -  1 

oooou 

VCOU  *  NCOi.  ♦  1 

000020 

1  F  <  3  •  !H'  35  ,  30  .  30 

000022 

30 

HAL 9  *  * 

000023 

h  •  hZ02 

000025 

CO  TO  43 

G  0  0  fl  25 

35 

H  G  A  H  i  -1 

00002ft 

40 

no  42  J«1 .  6 

0000 30 

JO  42  <1  *  1.2 

000031 

15  «  S  5  T ( <  l  )  ♦  40«*(J»1) 

07003ft 

l_E  *  V  r  IH  <L  )  *  40?*:j-li 

000044 

U0  42  1  «  LS, 16,10 

00004ft 

«U( ! ♦ 1 )  *  »D( 1*3) 

000050 

8Q< 1*2)  *  PD< 1*4) 

000051 

RO(I*3)  «  «D<!*5> 

C  0  0  0  5  3 

«u< ; *4)  «  net i*6i 

000054 

HU ( 1*5)  t  HQ( 1*7) 

0  1  C  0  5  4 

4? 

HO  (  1  )  *  P  O  U  ♦  9  ) 

0*0060 

R  9  T  o  9  *» 

c 

r 

Ks 

Sh!FT1vj3  OCCURS  4FTb3  0HtytnToo  C7rLu  j$  cf0HBlETc 

00006ft 

c 

4  • 

w*;_P  *  H*w9  .  i 

00087? 

f  T i  *A.9  -  2  )  55,  50  .  55 

000077 

>0 

H  •  HQ 

?  •?  0  '  7  4 

o:  TO  65 

000*74 

55 

T  *  T  ♦  h/ 02 

000076 

*0  TO  65 

000*77 

60 

T  *  T  ♦  H 

000101 

55 

n  ?  7  5  J  *  t  .  ft 

r  *  *  «  *  I 

0  **  7  5  <  L  *  1  .  ? 

030’?* 

{.4  *  v  5  T  <  <  l  7  *  4  6  ft-  •  {  j  - ;  ) 

0^0111 

w  *  «  %  i  O  (  <  L  ’  •  4  0  0 • f  j -  1  ) 

000117 

00  7  5  i  * l S  »  w fc , 1 0 

00012’ 

9(l*«)  •  t4  T  t  1*7) 

30012' 

-■'ll*?)  *  HO  <;♦{,) 

0  3  0  1  2  « 

-.'(I**)  *  h  ~  ;  i  *  5  > 

000126 

i  (1*3)  *  00(1*4) 

030127 

- : ( i  ♦  4 )  «  Hrii*') 

39913! 

80  ( S ♦ 3 5 

a 

SO  < 1 «  ? ) 

00013? 

80(1*2) 

a 

80(1^1) 

000134 

80(1*1) 

a 

brt(  I  ) 

900134 

M  80(1  ) 

a 

80(1*9) 

000144 

RETURN 

k  *  d 

C  THE  ffluOWINS  SUTEHENTS  EXIST  F OR  CH4NG I NG  ThF  STrP  SI2F 


C 


000146 

90 

iF(lXUk)  1 75 i 

135,  124 

009150 

124 

1*  •  -1 

000151 

125 

mPTN  a  nPTn  • 

l 

0  0  0153 

no 

I M  ■  M  ♦  1 

000155 

fNOP  a  o 

000154 

1:1  »  0,8 

000197 

VC3U  a  \iCOj  ♦ 

1 

000160 

return 

r  double  the  INTERVAL  if  possible 


000161 

000164 

C  0  0 1 6  6 

r 

135 

140 

!F(n  •  8 )  116 ,  140,  140 

K Z  «  «d*TN/2 

IF(2*K7  -  4PTN)  150,  125,  150 

c 

c 

I F  He'T'd  IS  ODU  TwgPg  IS  AN  EVEN  NU**BFa  Cf  PC’NT$  LFF*  to  CPNPyTfc 

000167 

c 

150 

1  r ;  (  <.  hnax  »  'M*,  ax  )  *55,  125,  125 

009172 

155 

IF  i N 0 T S  •  l)  125,  125,  160 

000175 

150 

NPTS  «  np T S/2 

000177 

mPTN  a  nRTn/P 

000200 

I**  »  5 

000231 

DO  165  J  •  1,6 

000202 

03  165  *L*1 »  2 

000203 

LS  •  NST ( <L )  ♦  400#(J*1) 

000210 

LE  •  N;D ( <L )  *  400*(J*1) 

000216 

DO  165  I  *  lS.lc.10 

000220 

80 < 1*1)  a  F o < 1*2) 

00022? 

80(1*2)  «  80 ( I *4) 

000223 

en< | *3)  <  HO ( 1*6) 

000225 

165 

BO(  I  *4)  i  80(1*6) 

050234 

H  2  ,  •  4 

000236 

* 

GO  TO  no 

r.  ULF  T^i  STEP  SIZE  IF  POSSIBLE 


r 


000236 

175 

t r  t  («d-miN)/MN!N> 

124, 

124,  189 

990241 

110 

N  •  al3j(E1/Enax 

*  1, 

)/. 6431471803 

00025? 

03  l 65  I  « 

1,N 

000252 

*  «  1 

000253 

I  r  <•»/?,  •«< 

*  -iH !  s »  195,  198 ,  105 

00026? 

195 

C3NT I M JE 

000765 

19? 

d  •  s 

000267 

*o  to  ?:o 

09026* 

196 

N  •  5  -1 

000271 

200 

*AlP  »  4 

005272 

a  J 

090273 

IF  <  I«  '  4  ) 

fl*. 

213, 

210 

885275 

71? 

T  a  T.  4 

000277 

9>iT  •  T 

00030" 

w  «-  1 

030501 

10  T?  21- 

090*9’ 

T  "  «  4  T 

090*94 

m3  T v  t  “  P  T  \  *  4 

0  3  0  3  C  5 

09030* 

?  1  4 

■J  9  0  3 1  * 

w07%iM7T9i*’*^4' 

: 1  ‘  120 

o  D  22?  J  *  x . 6 

(-"  0  32? 

no  ?2?  *L » 1 . 2 

090*23 

L S  •  \iT ( *L  )  * 

000530 

L  f  *  V;f:(<L)  ♦ 

G  9  C  3  3  * 

j:  ??2  i  *  os 

000340 

K  1  !  *  \\ 

0  g  C  3  42 

?2? 

t*( n  •  *r(K i 

090353 

•4D  *  h/?,*»N 

090157 

-/  2  s  *0/2. 

0*0361 

-  =  *20? 

09036? 

J  H  *  0 

000142 

isns  *  9 

oae363 

9fcTU?N 

000364 

END 

SUPROJHvt  :J* 

0  9  0  9  0  ? 

n  J  *  t  N  S  1  3  M  V  S  6 ' 

0  5  0  C  C  ? 


**09;? 
: :  c  5 1  ? 
oioco  * 

0509:4 

m  L’  >  3  k 

0  5  0  5  0  * 
0  9  0  9 O’ 
n  a  2 1*  i  ? 
990014 
0  9  0  51*- 
090-3' 
C'O'jf 

;  o :  o  * ' 

**  *  *  m  ty  * 

^  A  *  *  .»  C 

.  nt  *  J 

0  9. *  *  ? 
5  9  J  ?  s  ? 


1  fc-MIQ.  2C.  2),  V  ( 1 0  •  2  3  •  2  ) ,  -V  1 1 5  ,  2  0 , 2  ) ,  Ti(lS.2n,2)  PSC’0.2). 

2  v  D  (  2  9 , 2  >  ,  DNl<2),8P.t0,2".12).T.*fl9C05.rr(*0),M“(2'') 

3, :t(10.29.2)  i fe“ D” T (  ? ) 

4,  -h|NT{23;?) 

£  OU!  V*  _fcN:e  (TAMM.nSEM  ?  <TA(47<9)  ,£4j>.e)  (04(139). DO. 

1  (d9(40i>.DH) ,  c b r  1 1 > , fe - j ,  ?arHi6?'i).v).  (9*<?ooii,pyM, 

2  (TA(55).F*),  (  T  4  (  4  7a  )  #  (,),  (  r  r  (  3  > .  D  *  > »  <80(401),  TL  *  . 

3  ( 0i  (  95) .  ?S> .  (  H  <  31 1 ) .  (  T  4  (  52  ) «  Dvri,  (  r  c  (  j  > ,  *  i . 

4  (04(4*4).  ?*oG)  ,  (  U  <2’)  .G) .  (TA(3P),9«,(T4(?i»n,94T)J  (T4(9 

5  V  4 T  ■» )  ,  (  T  4 ( ? '  )  , 9  4  T - s  )  ,  (04(3  2 )  ,  U4T*x> 

6.  i  -0(1201)1  )T )  .  (  t4(  36)  .  f<rOT  ) 

7,  (04(54), ;“!'**) 

d'h-ov  /fc/^o.N**  /■o-‘D/c/*4 
SS-P(I)  *  9. 

c\-p  <  2 )  *  9 , 

S  v  “p  *  Jl 

<:  -o»  *  9, 

;  *  ?  0  0  !  «  1.2 
f  4 1 1 )  «  1 4  d  n  ’  t :  > 

1  «  v  5  ?  r  <  r ) 

. '  ? ;  o  j  *  i ,  > 

•MHH.J.M  .GT.  ,  ”  !  .  .  ’  )  )  'D  T 

•  i  *  c . : 

■ . ( i .  j ,  '■  >  *  :.c 

.(!._.!)  -  ' . 9 

*  -  (  1  : .  -  .  t  )  1  9  ,  * 


ooom 

uvuo#j«n  •  o,o 

685106 

a*<  i.j.n  ■  o,o 

cams 

DH?  1,J,M  *  0,0 

080120 

DT (  1*J,1)  •  0,0 

58312? 

NH  (  1  ♦  1 6 )  *  J*10*l*200*<  t •! ) 

080134 

ascj«n  *  o, 

5CG1S7 

V3<  Ji I <  «  0, 

U98t  43 

ontu>  »  twn>  •  tp<j,n 

230156 

ON C  J •  I )  «•  0, 

005154 

20 

CONTINUE 

080154 

C-NMP(I)  ■  SMHPUJ  *  DMJ.n  *  nH(UJ.I) 

C80U6 

SUMP  «  SUM"  ♦  DN(J.!>  *  BNfitJd)  * 

030207 

*  SU4PA  *  OV 1 1#  J»  ! )  ♦  ENIliJil)  *  DNtJ 

Of 022? 

FA  4  !  5  •  F  A  «  1  >  -  DNIJiP  »  EHdiJiP 

03023S 

r A  « l  >«A4AX1(FA( J ), 0,49E«10> 

•'68240 

200 

CONTINUE 

008244 

A  K  *  A  {  X  > 

058246 

A*  =  A*  -  AP<XJ  *  OX 

mi5S 

A2  •  AX 

083253 

AJOA2  ■  A1/A2 

C0G255 

PrtIP  *  -<EVH4{1)  ♦  ENMP<2M 

050257 

PH!  •  s' A  ( 1  >  *  F  A  ( 2  i 

033261 

U  •  »H!  /  4H0G  /  AX 

03026* 

PI  «  U  *  (*P«IP)  *  D*  /  ' Q» A2 7 

533771 

P?  ■  «SJPP  /  !Q*A2> 

038273 

P3  ■  «5 JPP *  /  (G«A2) 

060275 

P4  ■  *R40G  «  U**2  •  ( A 1/ A  2  -1,0)/Q 

OOOSOS 

p  a  »1»32*»3*P4  *p 

000310 

»AT«EM3DT(i:/EhD0T(2) 

ooc:i2 

IF « FA  <  n*FA<2)  ,  QT  (  1,E*10>  RATpFA  f  1 ) /FA  ( 2  > 

00031V 

RAlP  ■  (FA<J,)/FA(2)*ENMPf2).ENMPm )  t  FA(2) 

000S23 

! F  (RAT  ,QE,  3ATHN)  GO  TO  1262 

055329 

rat  ■  rathn 

098326 

RATP  a  3, 

030327 

1262 

IF  (RAT  .LE,  «AT*X>  GO  TO  1264 

Cf3S32 

rat  ■  RATHX 

C3S33? 

RATP  «  0, 

C053.J3 

1254 

CONTINUE 

030333 

RETURN 

0  0  0  S  3  4 

ENO 

74 


S  /B=»Q  JT  I  N=  TuPLCJiD 

010  305  ni^fc'iS?  JN  <K<2U#2)»Kt>lSC<2f'#2>.XfNTRI20»?)*WPNYR«20#?i«*8H,3»2>* 

1RS<20.?).T.I1&,20.2>,TWB!NT<20.2>.  DFIX  ( 20 . 2) .  v  ( 10 . 20 »  ? )  * 

2  t«C19.20.?>,DUM!TU?G.2>«nUMTU?0,2).nNUO,?0.2».»ETA<20.2>. 
3rNJM«23»2).DYUll*20,2),rC(*0>»W*«<20).BnC10t20,12>*TAlinoO) 

000005  FUJlVALiNC:  I  T A  ( 1?6  J , KK) ,  (  T A <  270  ) , * 0 1  SR i  •  <  Ta  ( 559) ,  XFNYM) » 

l(TAC72?>«w  =  NT9),  (TA(21j>»bW>.  I TAI96),RS>*<RD<801>.TL> 

1,  (TM759J.TWBINT),  <  T  A  ( *  Q9 )  *  DEL  X  ) » I »0<  16« 1 ) ,  V  > .  <  RP  ( 1  >  #  EM ) . 

3  (  T  A  <  84  9  > ,  OU M  PL )  *  (  Y  A  <  8*9  ) ,  0  J*  Yi  ) ,  ( 8DC  40 1 ) » DM  > »  t  T *  ( 6  47  > » RETA  ) , 

4  (Ta(433)»  •  NHM )  i  (PO(1?01),OY),  <FC(1)»X).  ( TA  ( 26 )  •  T  wO®  I  » • 

5  <  T A  <  4  7  5  >  *  CP*'.  <YA«474i,r<*V),  (rC(3)*RX) 


000005 

000005 

163 

1  !  *  <<<  Ji  P 

000011 

4  •>  |  SC  C  J»  1  >  8  1 

000015 

GO  TO  <164, 165 

00002? 

164 

<<u,n  8  l 

000026 

I F  (  X  )  601.  61)  1 »  600 

00003* 

600 

XfcVTRI  J, I )  *  X 

000034 

WF'JTRl  J,  1)  =  W W <  J,  1  ) 

000343 

XL  ■  X < ( «S < J,  1  > ) 

000052 

W  8  XL  *  TJOPI  *  RS(J.I)  •  RMO < TL < 1 » J* ! ) •  I  > 

000070 

Tl  ■  W  *  C 3 A  /TKAV  *  HTIRSIJ.P)  72,  /TRAP  1 7RS I J # I > 7RS < J . P 

03010* 

GO  TO  602 

000106 

601 

XfcNTRIJ.P  «  -.001 

03011? 

XL  s  X4(RS< J. 1 n 

080121 

wfcNTR(J.I)  *  XL*TwnpI*RS< J, !)*BHO(TL(l» J. P. ! > 

000142 

71  «  W  =  NTK<J,p  *  CPA  *  RT<RS< J. ! ) > /PKAV/2 . /TwOP ! /RS ( Ji i 1 /»S C J* I ) 

00016? 

602 

T*HIMT<J»l)  t  «17, 

000167 

HELXU.l)  8V<l.J,p*EM<l,J,p*CPL(TL<l.J.P.P*  CT«B  !  RT  ( J»  I  >- 

L  T  L  ( 1  *  J  *  I ))  /TnHP  J  /RS(  J#  !)/P(40(Tl?1»J*  I  > .  t »  7562 , /XL 

000242 

0U"ITIU.I>  *  TLd.J.J) 

000252 

O  U  5  T  L  CJ*i>  *  YL  <  1#  J.  I ) 

000262 

155 

XL  «  X  < ( PS  <  J i 1 ) ) 

000271 

u  *  XL  *  T  40P I  *  RSI  J#  1  )*RHO(U  n.JiDit) 

000307 

77  «  M  *  C3A  'FKAV  *  HTIRSIJ.P)  7?,  7TW0P  l  7RS  I  J*  I  >  7RS IJ,  !  > 

000325 

T»  R  »  617, 

000326 

X'RA  t  XfeviTRI  J,  !  )  *  DfcLXI  J, ! > 

090335 

fF(X-XT9A)166#167,l67 

000340 

166 

COMTIMJs 

000340 

W  *  4fc NTR ( J » 1  )  ♦  <w.*fcNTR( J, |> >  /  ngLX(J.I)  *  IX-XFNTRf J. P > 

000360 

T»4  1*  OjHITLIJiI)  ♦  ITWR-DUNITLI  J.PWPELX*  J.P*<X*XENTRIJ,p> 

000377 

n r c i .  j .  i  >  *  itjh-oumtli  j.  i  n/cx 

000411 

PUmTl(J.I>  «  TWB 

000415 

DMll.J.p  S  •W/VIlfJ.l) 

000431 

R&TAIJ.P  *  W/RSI  J.  P/FNUMI  J.P 

000*44 

UriUil)  *  4 

000450 

GO  TO  156 

000450 

15/ 

DKluiP  *  -»'VI«»JiP 

000464 

pril.j.p  *  IT4B-0UHTLI  J.  m/D* 

000475 

ihht.i  j. t  >  »  Y4B 

000501 

Rfe  T  A  1  J  i  P  3  W/RSI  Ji  P/PNU4I  J,  p 

000514 

156 

RETURN 

030515 

PNt' 

75 


006009 
,  000009 


£ 

r 


000909 


SJ8R3UIINE  VAPOR(J,D 
COMMON  /l/rii9M/8/»0/C/T A 

0 1  MENS  I 3N  FBIOO),  NM(2U  S  ,  10(10110*12)  *  YA ( 1000  ) , MC3N( 2) 

1  «PE(2i «  B6U(20,?',  RS(20»*>.  FNUN(2e,2)»N<10, 20,21, 0Y(10,?8#2) 

2,FNyM( ?0,2) * ?L<10«  20,2),  EMUO»tO*f).  VM0«20,2) 

3  *09(10,29,21,  MW( 20,2) 

EQUIVALENCE  (TA(25?>,  RAT),  (TA(30),  P),  < TA (54 > , WCON) 

1, (TA(697)»  Pfi),  <TA(473)«  OFOT),  (TA(476)  ,PA) ,  (TA(047)  ,BEtA) 

2, (TA(  90),  R8 ) »  (TA<433),  FNUM )  ,  (80(401), DM),  (80(1201),  DT) 

3,  ( TA ( 26 S  #  TN9P1),  (TA(393),7NUH) ,  (TA(474),  FKAV),  (YAI264),  T8) 

4, (80(801),  TU,  (80(1),  EMI, (80(1601),  V).  ( TA ( 475 ) , CPA  ) 

5,  ( Y A ( 492 ) ,  GUR),  <7A(J02),RV) 

6*  (TA(216),UN),  (80(401),  09) 


000005 

•E(l)  *  0,0 

5 

i 

i 

00000 A 

PE<2)  ■  0,0 

000007 

BET  A ( 1 1  »  WCON(!)*OR8T*ALOG<(P-PE(T>)/(P«PA)> 

000029 

M  «  BETA1 J, 1 )*RS{ J, I)*FNUM( J, I) 

? 

000040 

IF  (W)  190,  192,  192 

! 

000041 

190 

M  ■  0,0 

l 

00004? 

MHITE  (6,  2050) 

; 

000046 

2050 

FORMAT (  19MPW  NgG ,  SET  TO  ZERO) 

f 

i 

090046 

OM(l,J,M  ■  0,0 

030054 

DTU.J.Il  ■  TWOP!#FKAV*RS(J, !)*FNUM< J, I)/EM(t, J, !)*(YG  -Tl ( 1. 

;* 

ij»  in  /  c»l<tl<i,  j«  i).  f )  /  v(i,  j,  n 

I 

009125 

MM ( J,  I)  *  M 

| 

* 

000131 

RETJRN 

900132 

19? 

OMll,  j,  I)  *  -M  /  V(l,  Q,  |) 

1 

000146 

MM ( J t  I)  *  M 

i 

00015? 

Z  •  W/TMOP!*CPA/RS(J,II/FNUH(J,n/FKAV 

000164 

EX  ■  =XP  (Z)*l, 

000170 

niF«  AlM(Tl(1,J#I  ),I)  *  CPA*( T6*  TL ( 1 . J, I ) .  'FX  ♦  G JR* ( M/RV/RS < J,  l ) 
1  /RS(J.I)>**2 

000220 

CPLTL  •  CP^(TL(1,  J,  I),  I) 

000230 

200 

FORMAT  (1H1) 

000230 

OTd.J.H  »  •OlF/CPL(TLd.J,I),I>/EM(l,J,I)*W/V(l,J,l) 

000264 

IF(DT(1,J,I)  ♦  100,)  7711,  156.  156 

000273 

7711 

DTd.J.I)  •  -100, 

000301 

150 

RETURN 

T 

00039? 

END 

it 


] 


i  a 


i 

? 

I 


4 

i 

$ 


r 

c 

000003 

000003 


000003 

000004 

09000A 

000012 

000014 

000014 


function  x<(»s) 

CALCULATES  dissociation  rurnino  RATE  CONSTANT  as  FUNCTION  of 
RADIUS 

MANSION  XXK ( 3  ) ,  RRS  ( 5 ) 

data  XX 4 ,  RR8/  3.2,  6,3,  7.3,  0,2,  11,.  0  ,,  100  ,.  200  .,  300  .,  500  , 
/ 

CONVERT  RS(K)  TO  MICRONS 
RST  »  IS 

RST  ■  (ST  *  2,9 4E*4 

CALI  LN3RIN  (RST.XK.RRS.XXK.DXXDRS.S) 

XK  ■  X<  *1,6-3/2,54/2,34 
RETURN 

FNO  _ 


I 


4 


00000? 

000003 


SU'HOUUNfc  Pi.!'TY 

RETURN 

FNf 


000007 

000007 

000007 

000010 

000013 

Q10CH 

000023 

000024 

00002$ 

000055 

000054 

000073 

000074 

000075 

000077 

000101 

000104 

000110 

000114 

000115 


StHmuTJNfc  COMt>  ( oFU ,  P,  y,  MSP) 

COMPOSITION  OF  COMBUSTION  IU$fc$  AT  AniAQiTIC  FLAME  T^stPATuPg 
AS  A  F  jyrT  T  O'l  OF  O/F  PAT  1 0 ,  MOEBSilRP 

COMMON  /CM3T AH/  VSPtiS*  15,  3)  /PRSSP/  PBSSI3)  /OF /  M15U  ■<* 
fllMfcNSION  Y(NS®)«  VSPB(15,  15) 

VS  ■  NS 9 

IF  (0  ,.E,  PRSS  ( 2 )  )  INI)  «  1 
IF  (3  , 3T  *  PBSS525)  IND  t  ? 

A  e  <P  -  PSSSIINH))  /  <PP SSUNf)  *  1)  -  PRSS{fNf5)> 

HO  1  J  *  1,  N« 
no  1  I  *  l»  NS 

1  VS°«U,  1)  *  YSPIJ,  !,  IND)  *  A  •  f  YSP  ( J*  I,  IND  ♦!  )»YSP«  J,  I « JNO) ) 

no  2  1  •  1.  NS 

7  CALL  LNS^lN  <OFH,  Y(I),  Rf  YSPB«1»  I).  DYOX,  N») 

SUM  *  o. 

DO  -5  l  *  1*  NS 

IF  ( Y ( I )  ,36,  0,)  uQ  TO  S 

Y  (  l  )  •  0. 

3  SUM  *  3JM  *  Y ( 1 ) 
no  4  I  *  1,  NS 

4  Y  (  I  )  «  YU)  /  SUM 

RETURN 

pND 


000005 
000005 
00301 Y 
000021 


function  T3l<p,I) 

«o!L!N3  TfcMP  nr  Ox ( I  *  i),  rugl ( !  *  ®) 

r3"M0N  /PVf/  PV  { 15,  ?)  /TLX/  TL  <15.  2)  /QTl/  ''T|.(?) 

fALL  LNjRl  n  £  s  •  T«L*  PVIl,  I),  Tl?1,  I),  DTBDbW,  NTH!)) 

RETURN 

CNjl! 


093005 

000005 

300017 

000021 

000021 


function  viscvu.i) 

V6BOR  n  I  SC  3S I  T  y  FOP  nx»j  *  l),  FUEL  ( I  *  ?) 

COMMON  /TF.MX/  TFILM(15,  2)  /VlSCV*/  VV!S(15,  ?)  /DIPl*/  NTFLM ( 2 ) 
CALL  LN3'UN(T,  VISCV,  TFIl*(1.  D,  WISH,  |),  nVfSOT,  NTFlM<I>) 

vfscv  *  visrv*i,ofc-Y 

RETURN 


i»C  V  . 


FUNCTION  AMO  ( T|  I ) 

C  LIQUID  OENSITY  OF  0*U  ■  U,  FUEL!  t  ■  2) 

0 0 0 0 e !»  COMMON  /TL*/  TL <  Hi  2)  /RH6Y/  RH0LU5*  2)  /DTI/  NTL <  2) 

09000!!  SAIL  L M3W I M < T «  RMO»  TL<li  fl*  RMOLU.  U.  ORMODT,  NTUDI 

00001!  RETURN 

000021  END 


FUNCTION  AVAR  (T,  I ) 

C  VAPOR  *AESSU»E  OF  OX C  |  ■  1),  FUELU  ■  2) 

CCQ609  COMMON  /TL*/  TL(15i  2)  /AVY/  AVItS,  2>  /DTL/  NTL (21 

000005  CALL  LN3*lMUi  PVaR,  TL<1.  !).  PVU,  I),  DAVDTL*  NTL  U  > ) 


000017 

030021 

return 

END 

c 

FUNCTION  F<A  <T,U 

vapor  thermal  conductivity  of  o*u  ■  n,  fuelu  *  ?> 

000009 

REAL  KVAP 

000)105 

COMMON  /TFlHX/  TFILMU5,  2)  /KVAPY/  KVAPU5, 

2)  /DTFLM/  NTFLM(?) 

000005 

CALL  LN3RIN(T.  FKAi  TFIL*<1.  P.  KVAPU  I), 

DKVDTF,  NTFLMUI) 

000017 

rKA«FKA*l,0E-06 

000021 

return 

000021 

END 

000005 
000009 
000*1 » 
000021 


function  cva*  <T,n 

VAPOR  SPECIFIC  OF  OKU  ■  u,  FUELU  »  2) 

COMMON  /TFtMK/  TF  ILM(15*  2)  /CAVAf/  CAVAA ( 19,  2)  /nTFLM/  NTFLM<2) 
CALL  LNjRI N(  T«  CWA*,  TFJLMIl,  !),  CAVAA(  1#  D,  DCVDT,  NTFL*MI>> 
RETURN 
END 


000009 

000009 

00001Y 

000021 


FUND . jH  C»L  IT,  I) 

LIQUID  SPECIF t C  meat  OF  OXU  t  l), 
COMMON  /TLX/  Tl ( 1 5 i  2)  /CAV/  CP(15 
CALL  LN3AIS(Y»  C?L»  TLU*  U#  CA(1 

return 

FND 


FUELU  •  ?) 

2)  /DTL/  NTL(2> 

n.  dcadt,  ntlid) 


000005 

000009 

OQQOQ9 

000017 

000021 


78 


FUNCTION  Aj.M  U, I) 

MEAT  OF  VA»0AIZATJ0N  OF  OX ( I  *  l),  FUEL U  ■  ?) 

REAL  LAI 

COMMON  /TLX/  Tl ( 15i  2)  /LAMY/  L  AM ( 19 «  ?)  /DTL/  NTL  <  2 ) 
EcThu-'*3*1  Hi  T  '  TLUi  ! ) »  L  AM  ( 1 ,  I),  DALMOT,  NTLUU 

END 


09060? 
09000? 
30000? 
000016 
090029 
G  909  91 
00063!! 


00050? 

00000? 

00000? 

ooone 

000029 

000031 

00003? 


SUBROUTINE  NTM’NT  «OFH,P,KM,nOF RWM> 

MOkfecU'. *1  or  combustion  sases  as  a  function  or  o/r  ratio 

COMMON  /NTM/  NMTAK1I,  3)  /PRSSR/  P&SH3)  /or/  RI19),  NR 
OIM6NSUN  <(M T ( 3 J 
no  1  i  Si,  3 

1  CALL  LV3«H(OfR,  WMTU>«  R#  MMtAB(l,  | )  *  DUNMOF,  N») 

CALL  LN  j*(I  VI  P,  NM#  pass,  NwT,  dofrwm,  S) 

9ETURN 

PNO 


SURR3UT  I  Nfc  TGINT  < OFR , PC, ? , DTOOFR ) 

ADIABATIC  rL AME  TEMPERATURE  as  A  FUNCTION  OF  O/F  RATIO  ,  POESS, 
COMMON  /TGI/  TCK19,  3*  /PRSSR/  PRf5(3>  /OF/  R(19),  NR 
DIMENSION  TT { 3 ) 

DO  1  I  ■  1,  3 

1  CALL  LN3R I  N  (  'FRi  TT  1 1)  i  R,  TClfi,  I).  DTCDOF.  NR ) 

CALL  LVSRIN'PC.  T,  PRSS,  TT,  DT^OFR,  3) 

RCT'JRN 

END 


79 


000019 

cem? 

ooooi? 

000013 

000015 

000035 

00003? 

000040 

000045 

000051 

000053 

000055 

00900? 

00000? 

00011? 

000120 

000132 

000130 

000141 

000142 

090144 

000104 

000109 

000104 

oooio? 

000174 

00020? 

000204 

00020? 

00020? 


SUBROUTINE  VftSQAS  (OFR,  P»  T,  V|SQ(  TKR*  NS?) 

C  VI3C3SITV,  THERMAL  CONDUCTIVITY.  AMD  SPECIFIC  HEAT  OF 

e  COMBUSTION  OASES  AS  A  FUNCTION  OF  0/F  RATIO*  PRESSURE*  *NP 

C  TiMRIRA'URS 

COMMON  /TA1T/  TTABUS)  /VKSGS/  VSP(15*  15).  CPTARI15.  15)  /BT/  NT 
1/N4/  MM(15) 

DIMENSION  SUH(15),  V  ( 1 9  >  •  V(15).  VR(15>*  WMRtiS),  PHI  (15)  *  YPmKIS 
1).  C*(13) 

NS  ■  NS* 

DO  1  I  •  It  NS 

1  CALL  LNjRlNtT,  Vt|>,  TTAR,  VSP(1.  !)»  DYDX*  NT) 

CALL  C3MPOFR.  P.  V,  NS) 

SUMVIS  •  0. 

DO  20  I  •  1.  NS 
*r  (YU)  ,lT,  .001)  80  TO  ?0 
f  SUM( I )  ■  o. 

DO  10  J  •  1*  NS 

YR  (J)  ■  ?<!)/  V(J) 

MHR  (J)  ■  MM(J)  /  HM(!) 

PHI ( J)  ■  (1,  ♦  SORT  (VR(J)a  SORT  (HMR( J) ) ) >**?/(?, RJR  *  SORT  ( 

1  1,  •  1,  /  WMR(UU) 

Y*H1 ( J)  ■  YU)  *  PHKJI/YU) 

10  SUMU)  »  SUH  U)  «  YPMI(J) 

SUMVIS  ■  SUMVIS  *  VU)  /  SUMU) 

20  CONTINUE 

vise  ■  sjmvis  /i«E6 

DO  25  I  ■  1*  NS 

25  CALL  LN3RIN(T.  CP(I).  TT A8.  CPTABU,  I).  DCPDT*  NT) 

SI3A  •  0, 

SI3tfa  0, 

DO  30  t  >  1.  NS 

SI3A  a  SI6A  ♦  Y(l)  •  CPU) 

30  SUB  a  SI  OR  ♦  YU)  •  MM(!) 

CP3  a  S IGA  /  SI  SB 

TK8  a  V 1 SG  •  <  CPG  ♦  2,404/  91GB) 

RETURN 

END 


80 


$u«n3umt  Dirru<oF«,p,T,vPO?.  r,nrii,»<SPEC) 


00001? 

ooooi? 

00001? 

00001? 

ooooi? 

ooooi? 


ooooi? 


00001? 


000013 

OOOOI* 

090016 

000017 

000020 

000023 

000025 

000027 

000031 

00003? 

000042 

000063 
00007* 
000075 
0000  76 
000104 
000106 
000111 
000116 
00012? 
000125 
000131 
00014? 
000144 
00015" 
00019? 
000152 


r.  this  sjbpqutine  calculates  thf  nirrusiviTv  or  o*  thru 

C  41«U«1)  AMD  FUEL  THRU  M t X  U«2> 

c 

COMHON/:/  TA (1000) 

DIMfcMSfOM  4TM0u(?}.  0HPQT(?9)#X<T0PT<?9)*y{15) 

puiival=nc=  (  *tmol<i>,tam8> > 

CO'*H3N/Oir«,uS/EPSMl5,?>,8taM4U5,?),|PSOr(2)lS!80r(»)/Hl4/-H<15) 

DATA  AS/O/.npTS/"?/ 


nATA  0*:GT/ 

2,662  . 

2,318. 

? .  066, 

1,877, 

1.7?6, 

1,612. 

1,517, 

1,439. 

1.32, 

1.233, 

1,167, 

1,116, 

1,075. 

0,9996, 

0 , 949 , 

0,91?, 

0,66.16, 

0,961, 

0,8422, 

0,61?4, 

0.7712, 

0.7424, 

0.*640, 

0.6232, 

0,596, 

0,5756, 

"  , 5596, 

0.9352, 

0,517/ 

DATA  X<T06T/ 

.3. 

.4, 

.5. 

,6, 

,7. 

.8, 

1,0, 

1,2. 

1.4, 

1,6, 

1,9. 

2,0, 

2,5. 

3,0, 

3.5, 

4.0, 

9,5. 

5,0, 

6,6, 

8.0, 

10.0, 

20,0. 

30,0. 

irtKSiioo.io, 

40,0, 

100 

50.0, 

60., 

60,0, 

100,0/ 

r  cal:  const anst  for  this  propellant  SYSTPm 

c 

io  ns«nsp?: 

KS«1 
'•S*VE*  l 

no  20  <«i.ns 

PPSI(K.2)*=PS1(K«11 

Sr.H*(<,2)«StGMA(K,l) 

20  CONTINUE 

DO  50  J«1.2 
no  50  <»1 • MS 

PPSKK,  J>*  SORT(PPSorU)*EPSI(T,j)) 

SIGH*  (  <1  J)*(  o  |5*CStG0F«  JUSIGHA  ( K,  J>  ))**?/ 0,00422168/ 
1  SORT  (  t,0/MTHOL(J)  ♦  1.0/4HIKM 

50  CO^TINJ; 

100  CALL  C34PnF»,®#v,NS) 

SU“«D*Q , 3 

YN  JH«1  ,  3  *  VP02/P 
Tl*T/l, J 
T32*T1»SQNT(U) 

no  300  <«i.ns 
V|4)*V(<)*(P-VP0?>/P 
P«  T1/?»SI(«,I) 

r. ALL  LM3PlM(P.OHEG«XKTOFT,PHPGT»DYnx,MPT^) 

»!t*T32/3*<63/SlGHA(«»,  |  ) 

C.lNDtSjMU  ♦Y(4>/PD 
300  CC^INJ: 

nrj*  imuh/$u*D/p 

Pfc  To 

PNn 


8  i 
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Mixing  and  separation  criteria  of  hyporgolic  propellants  are  characterized 
by  a  chemical  reaction  time  and  a  mixing  time.  Separation  limit  is  defined  by 
assuming  the  two  characteristic  times  to  be  equal.  Experimental  results  with 
hydrazine  and  nitrogen  tetroxlde  propellants  agree  very  well  with  this  separa¬ 
tion  limit,  and  also  show  that  the  effect  of  chamber  pressure  Is  not  significant. 
MMH  and  A~50  with  rltiogen  tetroxlde  give  separation  limits  at  higher  pro¬ 
pellant  temperatures.  It  appears  that  stream  separation  always  occurs  by 
Impinging  compound  A  and  hyd'uzine  with  oxidize,  and  fuel  temperatures 
varying  between  C^P  to  30°F  and  40°  F  to  60*^’. 

A  second  work  area  of  the  project  was  concerned  with  the  extension 
of  the  Steody-Strte  Combustion  Computer  Mode*.  The  computer  program  w-r  s 
generalized  to  allow  increast-d  ver*  «v.iiity  in  running  various  propellant  systems. 
The  numerical  procedures  and  organizations  of  the  program  were  substantially 
lxnpro>  _d.  Propellant  property  values  were  gathered  for  several  systems  of 
Current  Interest.  A  computer  program  listing  is  presented. 
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